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1. 


Introdjction 


The  missile  map-matching  problem  for  guidance  updating  or  target  homing  is  shown  in  Figure  1,  The 
problem  as  defined  here  consists  of  locating  the  position  of  a  sensor  image  relative  to  a  reference  map 
which  is  stored  onboard  the  vehicle's  computer.  Once  the  match  location  is  found  the  rulative  location 
between  the  two  map  centers  can  be  used  to  update  the  vehicle's  navigational  position.  The  two  important 
\  performance  considerations  are  the  avoidance  of  false  fixes  as  measured  by  tiieir  frequency  of  occurrence 
V^^and  the  accuracy  with  which  the  position  fix  can  be  made. 


This  paper  describes  the  overall  design  and  evaluation  of  map-matching  systems,  (For  additional  back¬ 
ground  information,  the  reader  is  referred  to  references  1  -  10. )  Figure  Z  shows  th6  layout  of  a  typical 
system  design.  Here  the  system  parameters  associated  with  the  reference  data,  sensor  environment  and 
vehicle  are  integrated  to  determine  a  model  of  the  image  dynamics.  This  model  usep  in  conjunction  with  a 
signature  prediction  model  is  used  to  construct  the  reference  map  or  set  of  reference  maps  to  be  stored  on 
the  vehicle. 


In  the  map-matching  problem  a  number  of  errors  can  develop  between  the  sensed  image  onboard  the 
missile  platform  and  the  image  reference  map  stored  in  the  vehicle  computer.  Those  errors  can  be 
categorized  into  four  generic  classes  depending  on  their  impact  on  the  composition  of  the  sensed  image 
relative  to  the  reference  map.  Global  errors  which  impact  all  elements  in  the  sensed  image  are  generally 
accommodated  by  preprocessing  while  all  other  types  of  errors  must  be  accommodated  by  the  choice  of  the 
matching  algorithm.  The  scene  selection  process  is  important  for  determining  that  the  reference  map  area 
contains  sufficient  information  of  a  nonredundant  nature  to  successfully  perform  the  matching  task.  The 
scene  selection  process  consists  of  a  two  part  screening  process.  The  first  part  consists  of  various 
mathematical  tests  which  determine  to  a  first  level  the  amount  of  independent  information  and 
redundancy  within  the  scene.  The  second  part  consists  of  simulation  to  determine  the  acceptability  of  the 
scene  under  real  world  flight  conditions.  Finally,  3.  system  verification  process  is  required  to  determine 
from  the  nature  of  the  matching  data  whether  a  successful  match  has  taken  place  and  if  not,  what  appropriate 
action  should  be  taken. 

This  paper  is  divided  into  four  additional  section^^^^ection  2  describes  the  problems  associated  with 
describing  a  scene  mathematically  and  with  the  time  and  spatial-varying  nature  of  scene  signature  for 
various  sensor  types.  This  section  describes  the  nature  of  environmental  factors  on  image  dynamics  and 
their  impact  which  can  be  measured  in  terms  of  predictive  errors,  nonstructured  errors,  and  contrast 
reversals  for  various  sensor  wavelengths.  Finally,  remedies  are  discussed  which  can  mitigate  the  effects 
of  errors  due  to  image  dynamics. 


S^Section  3  describes  the  problems  associated  with  reference  map  construction  and  discusses  the  scene 
selection  process  by  which^^ood''^ reference  scenes  are  progressively  screened  out  from  those  that  are 
^ot  so  goo'a”, 

J 

‘^Section  4  describes  the  compatibility  of  various  classes  of  algorithms  to  accommodate  each  of  the  four 
categories  of  error  sources.  Ultimately,  since  the  magnitude  cf  these  errors  is  sensor  dependent,  this 
section  cros correlates  algorithm  appropriateness  for  each  sensor  wavelength. 

'^Finally,  Section  5  describes  die  mathematical  process  behind  developing  measures  for  system  perform¬ 
ance.  This  section  is  divided  into  two  parts.  Jn  toe  first  part  the  general  scene/error  model  is  discussed 
and  a  mathematical  approach  (through  a  list  assumptions  and  approximations)  is  outlined  which  can  be 
used  to  predict  the  probability  of  false  mat.ch  occurrence  based  on  a  number  of  system  parameters.  In  the 
second  part  of  the  section,  one  is  concerned  witli  (given  a  particular  scene)  using  the  statistical  data  from 
the  map-matching  algorithm  to  estimate  the  system  performance  in  near  real-time 
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Figure  2,  Overall  map-matching  system  design 
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lm4Be  dynamics  and  Its  Impact  on  comparison  of  sensed  Image  and  reference  ma 


Indicated  In  the  figure,  the  sensor  Amage  Interaction  Is  Influenced  by  a  number  of  environmental  factors. 
These  factors,  combined  with  Inherent  time-varying  material  physical  and  electrical  properties,  produce  an 
oscillation  in  the  scene  signature.  Dynamic  changes  In  the  sensor  scene  when  compared  to  a  time¬ 
stationary  reference  scene  can  cause  significant  errors  to  exist  between  the  two  maps.  These  errors,  if 
unaccounted  for,  are  generally  a  major  cause  of  failure  In  map-matching  systems. 

It  is  the  purpose  of  this  section  of  the  paper  to; 

1.  describe  the  scene  composition, 

2.  discuss  the  Impact  of  environmental  and  Inherent  scene  factors  on  signature  dynamics, 

3.  discuss  and  quantify  the  nature  of  the  map  difference  errors  when  sensed  and  reference  map 
are  compared,  and 

4.  outline  remedies  for  accommodating  map  difference  errors  in  the  system. 

As  the  influence  of  the  environmental  and  inherent  scene  factors  is  wavelength  or  frequency  dependent, 
the  discussion  will  focus  on  the  most  common  active  and  passive  sensor  categories  (1.  e.  ,  optical/near  IR, 
middle  IR,  thermal  IR,  and  microwave). 

The  following  section  will  describe  the  reference  m.ap  selection  process  including  methods  for  choosing 
reference  maps  to  reduce  the  map  difference  problem.  The  subsequent  section  will  discuss  the  role  of 
various  types  of  algorithms  in  accommodating  map  difference  algorithm  and  other  types  of  system  errors. 
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Figure  3.  The  impact  of  image  dynamics  on  map  comparison 
2.  1  Composition  of  the  scene 

The  scene  is  tlie  most  complex  component  of  the  map-matching  problem  and  the  most  difficult  to  model. 
Scenes  can  be  described  in  the  visual  domain  (the  eyeball  process)  as  being  composed  of  a  set  of  features. 
Actual  sensor  data,  broken  down  by  resolution  elements,  are  described  by  a  set  of  intensity  values.  There 
are  regions  of  Intensity  values  in  Ae  scene  which  can  be  considered  aralogous  to  features  in  the  visual 
domain.  These  are  homogeneous  regions*  within  the  scene  which  can  be  considered  equivalent  to  features 
(because  a  feature  can  be  defined  by  a  single  homogeneous  region  or  set  of  homogeneous  regions).  From  a 
physical  standpoint,  homogeneous  regions  are  areas  in  which  the  signature  (reflectance  for  visual  and  radar, 
emitted  power  for  middle  and  thermal  IR,  and  altitude  for  terrain  contours)  is  expected  to  remain  fairly 
constant,  e. g. ,  a  grassy  field  in  which  all  the  elements  in  the  region  are  expected  to  have  the  same  mean 
value  but  not  necessarily  a  constant  value. 

There  exist  variations  in  the  intensity  level  within  a  ho’^  logeneous  region.  Neglecting  the  possibility  of 
sensor  noise,  diis  signatur;'  variation  can  be  -attributed  to  eiUier  scene  resolution  constraints  or  texture 
variation  within  the  region.  Scene  resolution  constraints  can  cause  a  perturbation  in  the  signature  of  the 

e  We  define  a  homogeneous  region  to  be  a  set  of  spatially  connected  pixels  or  elements  which  possess  die 
statistical  property  of  at  least  first-order  stationarity  (constant  mean  intensity  level  over  the  region)  and 
possibly  second-order  stationarity  (mean  and  variance  constant  and  autocorrelation  independent  of  position). 
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regioiu  For  IniUnce,  one  cu  coneider  ^e  greeey  field  not  to  be  tinlform  but  Insteed  to  heve  e  few  fellen 
tree  trunke  end  ihrubi  dieperied  within  it.  If  the  ground  resolution  of  the  sensor  is  of  die  seme  megnitude 
es  the  sise  of  the  shrubs  end  tree  turnks,  then  we  would  ejcpect  verletlons  in  the  intensity  of  the  gressy 
regions  due  tu  these  objecU.  lit  should  be  noted  diet  if  die  resolution  of  the  sensor  were  to  increese  to  the 
point  diet  dimensions  of  objects  within  the  gressy  field  covered  severel  sensor  resolution  elements,  then 
these  objects  would  be  considered  homogeneous  regions  in  dismselves.  In  our  tree  trunk  exemple,  furdier 
increese  in  sensor  resolution  would  result  eventuelly  in  the  moss  on  the  fellen  tree  trunk  becoming  e 
homogeneous  region.  .Obviously,  die  process  of  identifying  homogeneous  regions  could  continue  ed  Infinitum 
es  the  sensor  resolution  wes  increesed. 

Thus,  we  cen  furdier  cetegoriee  e  homogeneous  region  in  the  physicel  domein  by  die  number  of  resolu¬ 
tion  elements  conuining  objects  which  contribute  to  e  sigceture  verletlon  end  in  the  stetisticel  domein  by 
the  number  of  sutisticelly  indeperdent  elements  which  comprise  the  region.  The  "scene  resolution"  con¬ 
cept  (11)  provides  e  useful  fremework  for  enelysing  the  steHsticel  verietion  of  e  region**.  We  shell  define 
this  scene  resolution  es  the  retio  of  the  everege  of  the  nur.ber  of  tensor  resolution  cells  to  diet  required  to 
meke  up  the  ecuivelent  of  one  independent  element  in  the  imeged  mep.  As  discussed  ebove,  sensor  resolu¬ 
tion  constreints  ere  one  contributing  fector  to  "scene  resolution"- -the  od.er  being  texture. 

Texture,  ceused  by  physicel  end  electricel  meteriel  verietlons,  cen  exist  even  within  purely  homogeneous 
regions.  The  three  primery  sources  of  homogeneous  meteriel  texture  ere:  Uluminetor-terget-detector 
geometry,  which  Includes  slope  end  slope  esimuth;  directions!  reflectence  end  ebsorpteice  described  by 
electromegnetic  theory  (Fresnel's  equetions)  end  surfece  roughness  effects.  Texture  produced  by  these 
processes  cen  be  virtuelly  resolution  independent  in  comperison  to  those  observed  within  e  resolution  depen¬ 
dent  homogeneous  region  (i.  e. ,  see  previous  discussion).  A  more  deteiled  presentetion  of  homogeneous 
meteriel  texture  '.’i  given  In  the  eppendix. 

b,  2  Signature  dviiemics 

In  order  to  estimete  the  intensity  megnitude  end  oscilletions  thet  occur  in  sensor  imegery,  it  is  first 
necessery  to  understand  tlie  relevant  physicel  end  electricel  meteriel  properties  end  governing  atmospheric 
end  meteorological  per/imeters  present.  A  summery  of  the  governing  meteriel  properties  for  each  sensor 
region  is  given  in  Teblu  1.  Similarly,  relevant  atmospheric  end  meteorological  parameters  for  each 
spectral  region  ere  gh'en  in  Table  2. 

Contrast  reversals  ere  of  importance  to  &e  mission  plen::er  because  of  the  potentially  decorreleting 
effect  they  cen  heve  e.n  mep-metching  system  performance.  A  summery  of  the  relevant  parameters  in  each 
spectral  region  thet  cen  Induce  these  effects  is  given  in  Table  3.  The  diurnal  end  seasonal  impact  on 
reference  area  signdture  characteristics  is  also  important  since  it  provides  the  mission  planner  widi  a 
time-frame  estimete  of  when  region  level  shifts,  hence  contrast  reversals,  ere  likely  to  occur.  A  summery 
of  0»e  time-cycle  impact  on  reference  area  signature  characteristics  for  each  spectral  region  is  given  in 
Table  4. 

The  impact  of  physicel  end  electricel  meteriel  properties  end  atmospheric  end  meteorological  effects  on 
time-varying  reference  area  signature  characteristics  will  now  be  presented  for  each  sensor  region.  The 
Impact  of  snow /ice /water  on  die  reference  area  signature  will  not  be  considered  here.  An  estimate  of  die 
magnitude  of  contrast  reversals  it  can  induce  within  typical  reference  areas  is  given  in  Section  2.  3. 

Pas^iive  ontic^/near  IR.  The  governing  msterial  and  atmospheric  properties  in  the  passive  optical/ 
near  IR  interval  (.  4^  ->  1.  6^)  are  short  wavelength  reflectance,  incident  irradiance,  atmospheric  attenuation, 
and  path  radiance,  respectively.  Contrast  reversals  in  this  spectral  region  are  primarily  due  to  changes  in 
material  reflectance  due  to  seasonal  effects  from  the  vegetation  growth  cycle. 

The  atmospheric  effects,  particularly  attenuation  and  path  radiance,  govern  die  degree  of  observed  con¬ 
trast  for  a  given  imaged  scene.  The  effect  of  atmrapheric  attenuation  is  to  uniformly  reduce  die  received 
radiance  across  the  scene.  Path  radiance,  however,  introduces  additive  energy  into  the  imaged  resolution 
element  via  direct  or  indirect  atmospheric  scattering  that  originated  outside  of  it.  The  net  effect  of  these 
two  terms  is  to  lower  the  observed  scene  signal-to-noise  ratio  (SNR)  for  a  given  sensor.  They  are  usually 
the  limiting  operational  factors  in  this  wavelength  region.  Complicating  operational  performance  pre¬ 
dictions  in  this  Interval  is  die  fact  diat  the  values  of  most  of  the  governing  material  and  atmospheric 
properties  are  generally  a  strong  function  of  v  velength. 

Passive  middle  IP..  The  governing  material  properties  in  die  paesive  middle  IR  {3p—5p)  rtgion  are 
middle  IR  reflecUnce  (hence  emittance)  and  thermal  inertia  (thermal  conductivity  over  the  square  root  of 
thermal  diffusivity).  T).e  predominant  atmospheric  properties  are  attenuation,  and 

♦  Presuming,  of  course,  diat  the  signature  of  die  objects  was  different  from  the  grass  at  the  wavelength  of 
..the  sensor. 

**  The  scene  resolution  is  computed  by  determining  die  number  of  Independent  (Nj)  elements  In  the  image 
and  then  dividing  dils  quantity  into  N,  die  total  number  of  resolution  elements  in  the  scene  (i.  e. ,  N/Nj). 
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Tabu  1.  Govaralni  Fhyaical  and  BUctrtcal  Malarial  Propartlaa 


(Decreaaittg  Order  of  Importance) 

Sanaor  Ration 

0 

Typa 

Pbyalcal 

flistiissL* 

Optical/Naar  IR 

P 

Surface  Roughnaaa 
and  Imaging  Geometry 

.4ti->I.SK  Raflectanca 

A 

Surface  RooghMaa 
and  Imaging  Geometry 

.4p->l.S  H  Ranectanca 

Middle  IR 

P 

Thermal  Inertia 

Imaging  Gaoittatry 

Surface  Roughnaao 

3  Reflectance 

.4H-.1.SP  Abaorptanca 
3K-*SH  Emittanca 

A 

Surface  Roughneae 
and  Imaging  Geometry 

Si*—  5  i*  Raflectanca 

Thermal  IR 

P 

Thermal  Inertia 

Imaging  Geometry 

Surface  Roughnaaa 

.4  it- 1.61*  Abaorptanca 

Si*—  12i*  Emittanca 

A 

Surface  Roughnaaa 
and  Imaging  Geometry 

Si*  — IZH  Reflectance 

Microwave 

P 

Surface  Raughnaaa 

Imaging  Caomatry 

ThermM  Inertia 

Microwava  Reflectance 
and  Emittanca 
.4i*—1.6i*  Abf'rptanca 

A 

Surt^ica  Roughnaaa 
and  Imaging  Geometry 

Microwave  Reflectance 

*  A*  Activa  ayatetn 
P«  Paaaiva  ayatam 

**  Diractloaal  alaetrical  propartlaa  esiat  io  aach  eaaa  which  vary  with  aariaca  roughaaat  and  Imagind 
(illuminator  •  aurfaca  •  aanaor)  gaomatry. 


Sanaor  Raaton 
Optical /Naar  I 


Mlddla  IR 


Tharmal  IF 


Microwava 


Tabla  2.  Atmoapharlc  and  Mataorological  Impact  on  Sanaad  Imagery 

Type  Paramatrr  Impact  on  Inaaary* 

P  -  Small  to  atrong  for  path  radiance  and  attenuation 

A  -  Small  to  atrong  tor  pad>  radiance  and  attenuation 

P  -  Small  to  moderate  tor  pati>  radiance 

-  Small  to  atrong  for  attenuation 

-  Small  to  moderata  for  raradiation 

-  Small  to  modarata  for  latent  and  aanaibla  heat  tranafar 
Oapanding  on  wind  epaad,  pracipitable  water,  and 
atmoapharic  and  ground  tamparaturaa. 

A  -  Small  to  atrong  for  attanuation 

P  -  SmaL'  for  pad)  radianca 

-  Small  to  atrong  for  attenuation  and  raradiation  dapanding 
on  apaciea,  concantratlon,  diameter  and  tamparatura  of 
aaroBol  diitributlon  pitaent. 

•  Small  to  medarata  for  latent  and  aanaibla  heat  tranafar 
depending  om  wind  apeed,  pracipitable  water  content  and 
atmoapharlc  and  ground  tamparaturaa. 

A  •  Small  to  atrong  for  attenuation 

P  •  Small  for  attanuation  unlaaa  rain  U  praaant 

•  Small  to  modarata  for  oxygen  or  water  (abaorption  and) 
raradiation  depending  on  apaciaa  concentration  praaant. 

•  Small  to  modarata  for  latent  and  aanaibla  heat  tranafar 
(Impacta  molatuTO  availability,  hence  material  amittanca 
and  tharnnal  energy  balance). 

A  -  Small  for  attenuation  unlaaa  rain  la  praaant 


a-  Aaaumoa  operation  under  cloud  cover  with  no  precipitation 

-  Atmoapharic  attenuation  ia  dependent  on  tt>e  apaciaa,  concantratloo  and  diameter  of  aaroaol  diatributiona 
praaant  and  atmoapharlc  preaaure  (govarna  molecular  apaciaa  concentration). 


—6“ 


T«bU  3.  Sourcaa  of  1im(o  Controa:  Savaroal* 


Sanaor  Raaion 

Zna. 

Cauaa 

Optimal/Naar  IR 

P  and  A 

-  Optical /Hear  IR 

Vegetation  Ranactanca 

Middle  IR 

P 

-  Material  Thermal  Inertia 

-  Diurnal  Solar 

Irradianca  Cnmponant 

A 

•  Middle  IR 

Vegetation  Reflectance 

Thermal  IR 

P 

-  Thermal  Inertia 

A 

•  Thermal  IR 

Vegetation  Reflectance 

Microwave 

P 

-  Atmoapherlc  Raradlttion 

-  Thermal  Inertia 

A 

-  Mlciowave 

Vegetation  iXaflectance 

*SBow/Ice /Water  conr^las  rnn  produca  contraat  ravaraala  in  aach  imaging  region 


TaUa  4«  Diurnal  and  Saaaonal  ^nvironmantal  Impact  on  Sanaad  Imagery 


Sanaor  Raaion  Type  Tima  Cycle  Impact  on  Imanarv 


Optical/Near  IR 

P 

Diurnal 

•  Small  to  atrong  (dependa  on  apoctral  nnd 
abaolute  level  of  lUuRvlnation  Imagery  la 
obtained  under). 

Saaaonal 

.  Small  for  apeetral  irradiance  ebangaa 
(aun'a  declination  angle)  but  moderate  for 
illumination  level. 

•  Moderate  to  atrong  over  vegetation  cycle 

A 

Saaaonal 

•  Moderate  to  atrong  ever  vegetation  cycle 

Middle  IR 

P 

Diurnal 

Saaaonal 

-  Strong:  chort  and  middle  wavelength 

irradiance  drivea  thermal  inertia, 
and  direct  reflected  Middle  IR 
component. 

-  Moderate  for  apeetral  and  abeoluta  Irradiance 
level  (hence  thermal  inertia)  differencea  from 
aun'a  decllaatloa  angle. 

-  Snail  to  moderate  over  vegetation  cycle 

A 

Saaaonal 

-  Small  io  moderate  over  vegetation  cycle 

Vherntal  IR 

P 

Diurnal 

•  Sroeg:  abort  wavalangth  Irradiance  drivea 
thermal  inertia. 

Saaaonal 

•  Same  aa  paaalve  Middle  IR 

A 

Seaaonat 

•  Snail  to  moderate  over  vegetation  cycle 

Microwave 

P 

Diurnal 

-  Small  to  moderate  for  high  microwave  omit¬ 
tance  objects  (tliarmal  inertia  can  dominate). 

-  Small  for  low  microwave  omittance  objecta 
(aky  temperature  dominatea). 

Saaeonal 

-  Small  to  moderate  for  high  microwave  emit  < 
tance  objects  (sun's  declination  angle). 

•  Snail  for  low  microwave  omlttanca  objects 

-  Small  to  moderate  over  vegetation  cycle  de¬ 
pending  on  caitopy  and  soil  moisture  content. 

A 

Saaaonal 

-  Moderate  over  vegetation  cycle  depending  on 
backscatter  coefficient. 

rcndlmtion.  Th*  priRcipil  m«teorolotic«l  intaraction  paramatera  ara  latant  and  aanalbfa  haat  tranafar 
(avaporation  and  convactloa  raapactlvaly).  Tha  matarial  tharmal  inartia  la  ralatad  to  tha  not  rata  of  haat 
axchanga  at  tha  aurface  batwaan  tha  air/matarlal  tntarfaca.  Tharmal  Inartia  offaeta  drivan  by  tiia  abaorbad 
abort  wavalacgth  Incldant  irradianca  during  tha  daytlma  ganarally  pradomtnata  tha  amlttad  powar  componant 
at  night.  During  tha  day.  both  raflactad  aolar  middla  n  and  tharmal  inartia  componanta  contribute  to  tha 
obaervad  aignature.  Contraat  ravaraala  in  diia  apactral  region  are  primarily  related  to  matarial  tharmal 
inertia,  where  the  amallar  the  magnitude  of  thia  parameter,  the  greater  the  temperature  (hence  amieaiva 
power)  oacillation.  Contraat  ravaraala  can  alao  be  induced  in  title  apactral  region  whan  a  large  aolar  and 
atmoapheric  3p  5fi  flux  ia  preaant.  coupled  with  a  low  to  moderate  nrtatarial  aurface  temperature.  Hare, 
the  time-varying  nature  of  the  downwelling  flux  paaaaa  through  a  cycle  of  amall  to  large  to  amall 
coincident  with  the  aolar  aenith  angle.  If  the  3p  —  Sp  aolar  raflactad  componant  ia  larger  than  that  due  to 
material  emiaaion.  a  time-varying  contraat  ravaraal  can  raanlt. 

The  limiting  atmoapheric  caae  due  to  attenuation  occura  whan  an  image  in  obtained  through  an  atmoaphera 
with  a  moderate  to  large  diameter  aeronol  dlatributlon.  In  thin  apactral  interval  atmoapheric  attenuation 
from  aeroaola  uaually  pradominatea  over  reradlation  from  praciptabla  water  during  the  daytime,  but  at 
night  the  raverae  ia  poaaible. 

Senaible  and  latent  heat  tranafer  can  impact  the  imaged  apectral  aignature  in  thia  region  >y  altering  the 
ground  temperature.  Thia  in  turn  impacta  tiie  image  aignature,  particularly  at  night  whan  i\a  ground 
emiaaion  component  predominatea. 

aive  thermal  iR.  The  governing  material  propertiea  in  the  paaalva  titarmal  IR  (Bp  —  K'p  )  region  are 
abort  wavelength  reflectance  (typically  .4p  —  1.6p  ),  tiiermal  IR  amittance  and  titcrmal  inertia-  The  princi¬ 
pal  atmoapheric  propertiea  and  meteorological  interaction  parametere  are  identical  to  tiioae  i.i  the  middle 
IR  region. 

In  tiiia  apectral  region,  material  tharmal  inartia  la  tiia  aole  cauae  of  obaarvad  ground  aignature  contraat 
reveraala.  Since  a  negligible  amount  of  thermal  IR  energy  emitted  by  the  nun  penetrataa  tita  atmoaphere, 
abort  wavelength  aolar  irradiance  driven,  thermal  inertia  effecta  predominate  the  image  over  the  diurnal 
cycle. 

Atmoapheric  attenuation  in  a  dry,  cloud-free  atmoaphere  ia  amall  in  title  apactral  region.  A  aubataKial 
amount  of  reradiation  (h^tnce  image  contract  reduction)  can  occur,  however,  whan  a  humid,  warm  atmoa¬ 
phere  ia  preaant  due  to  increaaing  e.niaaive  powar  with  precipitable  water  a^  atmoapheric  temperature. 

Such  condltiona  will  often  form  the  cloud-free  limiting  caae  for  aenaor  operation  la  thia  apactral  raglou. 

Senaible  and  latent  heat  tranafer  becomaa  important  whan  a  aignificant  differanca  in  atmoapheric  and 
ground  tamperature  exiata, coupled  with  a  non- aero  wind  apeed  a^  relative  humidity.  Thaae  haat  tranafer 
componanta  can  produce  a  noticeable  aignature  oacillation  for  a  reference  area  imaged  under  widely  varying 
meteorological  condltiona.  Furtiiermore,  the  magnitude  of  theae  parametara  are  often  difficult  to  evaluate 
due  to  the  lack  of  the  neceaaary  ground  truth  data. 

Paaeive  microwave.  The  governing  material  propertiea  in  tha  paaaive  microwave  imaging  (,3  cm  to  3.0  or) 
region  are  paaaive  microwave  reflectance  and  therixial  inei'tis.  The  principal  atmoapheric  parameter  hare 
ia  tite  contribution  of  precipitable  water  to  the  aky  brightneaa  temperature.  Senaible  and  latent  heat  trana¬ 
fer  componenta  tend  to  have  little  impact  on  the  obaervad  aignaturaa  unlaaa  high  microwave  amittance 
materiala  predominate. 


Contraat  reveraala  are  cnly  poaaible  in  thia  apectral  interval  ia  two  caeaa.  Tha  firat  involvea  materiala 
witii  low  microwave  reflactancea.  Here,  ttie  material  microwave  emittai'.ca  (tlmea  ground  tenaporatore) 
component  predominatea  and  the  raaulting  onergy  balance,  hence  imagery,  bahavac  aimilarly  to  tiiat  in  the 
thermal  IR  region.  In  the  aacond  and  much  rarer  caae,  a  ravaraal  will  occur  whan  tha  aky  brightneaa 
temperature  ia  greater  than  the  material  temperature,  Thia  ia  generally  only  poaaibl*  under  cloud  cover 
condltiona  when  a  ^batantial  amount  of  prec^itable  water  axiata  along  with  a  low  to  moderate  ground 
temperature  (~273”R -» -'Z90®K),  Hare,  the  emitted  energy  from  the  precipitable  water  becomaa  graater 
tiian  that  from  tiia  rafarance  matarial.  Aa  a  conaaquanca,  materiala  with  a  high  microwava  raflactancr 
(i.  a. ,  matal  and  water)  can  nave  greater  apparent  brightneaa  temparaturea  than  thoaa  with  a  high  micro¬ 
wave  amittance  (i.  e. ,  aoll).  Thia  raaulta  in  a  ravaraal  over  tha  expected  caae  where  a  dry  atmoaphere  ia 
preaant. 

Coatra:?y  to  general  belief,  the  only  materiala  that  can  not  exhibit  tha  firat  type  of  contraat  ravaraal  dia- 
cuaaed  above  in  the  35  GHa  and  94  GHa  banda  ara  metal  and  water,  ainca  moat  nnatariala  poaaeaa  high 
microwave  emittancaa  in  theae  ragiona  at  annall  acan  anglan.  Aa  a  conaaquanca,  regional  error  ahifta  (and  in 
aome  caaea  contraat  reveraala)  can  raanlt  alnce  many  common  materiala  (i.  a.  ,  vegetation,  aoil,  concrete 
and  rock)  exhibit  microwave  emittanee,  hence  tiiarmal  inertia  dominated  time-varying  oacUlationo.  Pre¬ 
diction  of  vegetation  and  aoil  aignature  magnitudaa  and  thair  oacillationa  can  be  vary  difficult,  however, 
bacauae  of  the  Impact  of  moiatura  availability  on  microwava  matarial  amittance.  Lika  in  tite  infrared 
legiona,  addition^  inatabUlty  in  the  microwava  aignature  can  occur  due  to  atmoapheric  reradlation  affacta; 
particularly  for  metal  and  water  which  poaaeaa  low  and  moderate  microwava  amittaocea  raapactlvaly. 
Paaaive  microwave  aignature  variationa  for  theae  materiala  ara  generally  much  larger  than  "in  the  infrared 
for  aimilar  condltiona  which  produce  atmoapheric  raradiation. 


Acttv  »v»t<nr>t.  For  Imaging  laaart  and  radars  Ihs  j^evamlng  *natartal  slretrlcal  proparty  is  raftsctaocs 
(or  tho  backs cattsr  cosfflcisnt).  Atniosphsrlc  absorption  and  scattsring  (attsnuatlon)  Is  oftsn  ths  limiting 
snvironmsntal  factor  for  lassr  Intaglng  cystoma,  al^ongh  It  Is  usually  small  for  radars.  Thsso  eystoms 
ars  at  l0v'.st  dlrsrtty  Insansltlra  to  many  of  tbs  comploa  tlme-varylug  physical  atmosphsrlc  and  mstoorologi- 
cal  sffscts  dtat  Impact  passive  systsms  (U  s. .  tharmal  Inertia,  solar  Irradlanca.  and  latoat  and  sensible 
boat  transfer). 

An  additional  class  of  active  sensors  existed  that  use  ttta  spectral  transmitted  beam  la  a  phasa  modulated 
carrier  or  range-gated  form.  The  advantage  of  these  sensors  Is  that  they  can  be  relatively  Insensitive  to 
all  material  and  meteorological  properties  and  generally  are  only  limited  by  atmospheric  effects.  In  ths 
first  case  a  frequency  modulated  signal  is  placed  on  an  optical  laser  carrier  beam.  Very  accurate  target 
x:.uging  and  depth  information  is  possible  by  detecting  the  phase  front  distortions  of  tho  retnxnod  beam 
induced  by  ths  object. 

The  second  type  of  system  is  operated  in  a  ranging  form  by  measuring  ths  two-way  propagation  time  to 
the  ground  or  target  (down  and  forward-looking  respoctivsly).  A  common  down-lookl^  fotm  of  dtis  system 
is  a  radar  altimeter  veed  in  Terrain  Contour  Matching  (TERCOM).  A  widely  used  forward-looking  form  is 
the  laser  rangefinder  used  in  tactical  armored  vehicles. 

For  both  types  of  systems  two  governing  performance  factors  sxist.  First,  the  reflected  object  energy 
must  be  high  enough  to  produce  an  acceptable  SNR.  This  often  limits  the  operational  distancs  because  of 
beam  dispersion,  and  atmospheric  effects  (lasers  only).  The  second  luvolves  the  beam  pattern  itself.  If 
it  is  too  large  in  diameter  versus  the  imaged  object,  phase  information  becomes  ambiguous  for  the  first  type 
of  system.  For  the  altimeter  or  rangefinder  system  this  also  poses  a  problem  due  to  an  increasing 
uncertainty  in  knowing  die  object  that  produced  the  first  p£  strongest  return.  Their  princVpal  disadvantages 
Include  hardware  complexity  and  lank  of  maturity  (phast*  modulated  laser)  and  potential  inaccuracy 
(altimeters  and  rangefinders)  against  point  targats  due  to  reference  imaging  requirements  and  their  usual 
one-dimensionri  configuration.  When  targeting  conditions  permit  and  operation  essentially  invariant  to 
environmental  conditions  is  necessa.'y,  these  two  types  of  sensor  systems  should  be  strongly  coneldered. 

Atmospheric  and  meteoroloelcal  parameters.  A  summary  of  die  relevant  atmospheric  and  metaorologi- 
cal  parameters  on  sensed  terrain  imagery  is  given  in  Table  2.  Ker%  molecular  absorption  has  not  been 
directly  considered.  It  is  at  least  im^led,  however,  since  the  atmospheric  windows  utilised  for  remote 
sensing  exist  in  regions  where  these  effects  are  smalL  Molecular  absorption  band  characteristics  vary  with 
temperature  and  pressure  for  a  given  species.  Aerosol  absorption  and  scattering  are  less  specific,  since 
they  also  vary  with  the  diameter  distribution  present. 

From  Table  2  it  is  apparent  that  Mth  radiance  affects  caused  by  aerosol  water  decrease  notieeably 
beyond  dte  optical/uear  and  middle  w  regions.  This  is  a  result  of  the  aerosol  diameter  dlstributioas  typi¬ 
cally  present  and  the  small  amount  of  sdar  irradiaace  that  exists  in  the  thermal  IR  and  passive  microwave 
regions.  Rerodiation  becomes  increasingly  important  wid>  wavelength,  and  in  passive  microwave  imagery  it 
is  die  dominant  rain-free  relevant  atmospheric  paratxatsr.  Latent  and  sensible  heat  transfer  «re  the  pre¬ 
dominant  meteorological  remote  sensing  parameters.,  and  can  have  a  moderate  impact  on  the  resulting 
energy  balance  present  in  middle  and  diermal  IR  imaging  and  alter  the  resulting  emittances  of  some  passive 
microwave  materials  (particularly  soil). 

Time-cycle  impact.  Large  oscillations  avd  possibly  contrast  reversals  in  material  signatures  often  occur 
during  diurnal  and  seasonal  time-frames.  A  summary  of  the  relevant  phenomena  for  active  and  passive 
sensor  systems  is  given  in  Table  4.  Two  factors  are  evident  from  die  data  given.  First,  die  performance 
of  each  passive  ssnsor  system  can  be  altered  by  the  level  and  spectral  distribution  of  incident  solar  irradi- 
ance  in  ths  atmosphere  and  at  the.  ground  plane.  Second,  the  spectral  reflectance,  diermal  inertia,  and 
moisture  availability  associated  with  vegetation  growth  cycles  on  land  can  significantly  impact  imaged  signa¬ 
tures  in  every  spectral  band  on  a  eaasonal  basis.  Only  phase -modulated  or  range-gat^  lasers  and  radars 
appear  to  be  relatively  immune  to  this  problem  as  long  as  dsciduous  tress  are  absent. 

Two  items  have  been  omitted  from  Table  4  for  simplicity.  The  magnitude  and  type  of  atmosphsric  and 
meteorological  effecta  praaant  prior  to  and  at  tha  momant  of  imaging  are  repraaentad  by  a  joint  diurnal- 
seasonal  tima  cycle  probability  distribution  function.  Likewita,  tha  prasance  of  eaow/ica/w&ter  widiin  the 
reference  aree  can  also  bo  dascribed  by  another  joint  diurnal-seaeonal  probability  distribution  function. 
These  two  dUtributions  are  very  complex  (perhepe  preeendy  indeterminate)  and  only  moderately  correleted 
wldt  time.  Coosequendy,  at  baet  it  ia  only  poeaible  to  approximate  die  impact  of  theee  factore  on  spactral. 
timo-varylng  reference  area  elgnaturae. 

2.  3  Mao  dUfarence  arrore 

From  a  syatams  point  of  viow  ona  can  categoriae  all  the  map  differencae  as  affecting: 

1.  die  spatial  shape  of  homogeneous  regions, 

2.  the  relative  mean  intenaity  Isvela  between  homogeneous  regiona,  and 

3.  tha  absoluta  intenaity  laval  of  a  roglon. 


In  tht  vnrnaculnr  ft***  aflnctt  arc  commonly  rotorrod  to  aa  nonatructurad  arrora,  eon'raat  rovaraala, 
and  pradictiva  coding  arpora,  raapactivaly.  A  eomUinatlon  of  theae  factora  ara  ganarally  praaant  In  aanaor 
Imagexy  and  Induca  arrora  in  tha  map-matching  procaaa  dua  to  thnir  complax  nature. 

Nonatructurad  arrora  can  be  brokan  down  into  two  catagoriaa.  In  tha  firat  eaaa,  tha  impact  of  tha  partur- 
bationa  ia  pradictabla,  although  it  may  not  ba  poaeibla  or  daairabla  to  prepare  a  large  number  of  rai'erence 
acanec  for  compenaation.  Errora  in  thia  claaa  include  ahadowa,  which  can  lead  to  contraat  raveraala 
within  tha  affected  region.  Their  location  can  ba  ealculaU-d  given  tha  illuminntor-target-vahicle  geometry 
combined  with  the  terrain  topography.  Errora  of  tha  avcond  type  are  more  diffic  ult  to  predict,  hence  to 
produce  companeating  imagea.  Thaae  errora  include  terrain  araaa  obocured  by  clouda  and  anow /ice/water. 
Hare,  fta  joint  probability  apaca-time  error  diatribution  affecting  the  reference  area  (hence  aach  reiolution 
element}  ia  virtually  unknown. 

Tha  net  effect  of  changea  in  tha  atmoapheric.  meteorological  and  phyaical  and  electrical  material  pro- 
pertiaa  ia  to  produce  variationa  in  the  intenaity  level  of  one  'hom,3grneoua  region  relative  to  another.  If  the 
intenaity  level  ahifta  are  aavere,  contraat  raveraala  between  regiona  can  raault.  An  aatimate  of  the 
expected  range  of  contraat  ratio  reveraala  between  repreaentative  natural  materiala  ia  given  in  Table  5. 
Maximum  valuea  and  the  governing  parameter  are  given  in  two  caaea  for  each  apectral  region.  In  the  firat 
caae,  contraat  revereal  rangea  due  to  phyaical  atmoapheric  and  meteorological  parametera  are  given.  In 
the  aecond  caae,  thoee  produced  by  a  now /ice /water  are  p  reaentnd. 


Strictly  apaaking,  aigoature  variationa  cauaed  by  anow /ice /water  are  predictive  errora.  The  effect  of 
Aia  complex  ia  to  produce  random  apace  and  time-varying  aignature  boundariea,  hence  artifical  homo- 
geneoua  regiona,  within  tho  reference  area.  Aa  a  reault,  contrart  reveraala  can  occur  within  the  aenaor 
image  due  to  aignature  variations  between  homogeneoua  regiona  created  by  the  snow /ice /water  and  tnoae 
from  the  nominal,  underlying  material  signature.  Preprocessing  techniques  that  emphaeiae  homogeneoua 
regions  in  the  aenaor  scene  can  produce  cataatrophic  ml^>-matching  failures  when  snow /ice /water  are  likely 
to  exist. 
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For  the  paeelve  optlcel.  neer  St  end  ell  ectlve  cesea,  the  output  le  given  in  db  ohenge  in  meteriel 
reflectence.  For  the  peative  middle  end  tb<'.rmel  IR,  end  micruweve  cease,  reeults  ere  preeented  in  both 
wetta /cm^.eteredien  end  db  of  redience  chenge  between  regiona.  Reeulte  ere  eimilerly  given  in  the  enow/ 
ice/weter  ceeee  except  the  aigiieture  of  the  perturbing  atete  is  compered  directly  to  e  nominal  meteriel. 
Reeulte  for  the  peeaive  middle  end  diermel  Ut,  end  microwevo  ceeee  were  determined  with  the  aid  of  e 
aophleticeted  atmospheric  boundary  layer  model.  Contrast  reversal  rengea  were  not  computed  for  man¬ 
made  materials  because  of  the  complexities  introduced  by  geometry  and  internal  heating  (for  the  middle  end 
thermal  IR  cases). 

Contrast  reversals  produced  by  means  other  than  snow/ice/waier  will  first  be  examined.  From  Table  5, 
it  ie  clear  that  the  vegeution  cycle  can  produced  significant  contrast  revarsa!a  against  soil  (as  well  as  other 
material)  backgrounds  for  active  and  passive  optical/near  IR  and  active  middle  and  thermal  IR  imaging 
systems.  The  largest  reversals  in  the  passive  middle  and  thermal  IR  cases,  however,  are  generally  pro¬ 
duced  by  solar  irradlance  driven  thermal  inertia  differences  between  materials  present.  Contrast  reversals 
can  also  occur  in  the  passive  microwave  region  due  to  the  vegetation  cycle,  where  the  primary  contributing 
factors  are  soil  and  plant  moisture  availability.  Reversals  or  intensity  shifts  between  homogenous  material 
regions  are  primarily  produced  in  this  spectral  interval  by  moisture  availability  and  thermal  inertia  effects 
which  Impact  the  microwave  emittance  and  gi  ound  temperature  (hence  foe  emissive  power  ground  component) 
depending  on  the  cIlRiatlc  conditions  present. 

When  contrast  reversals  due  to  predictive  errors  from  snow /ice /water  are  examined,  it  U  clear  that  foe 
magnitudes  produced  by  this  cause  are  greater  than  those  from  foe  corresponding  non-snow /ice-water  (vege¬ 
tation  cycle  and  thermal  inertia)  cases  in  every  instance.  Although  foese  valuea  may  sarve  as  reasonable 
uppei'  bounds,  foe  mission  planner  should  be  aware  that  changes  In  foe  snow/lra/water  state  can  produce 
aubstantlal  signature  valrations  over  a  short  to  long  time  frame  due  to  foe  complicated  physical  and 
electrical  properties  of  this  material  complex. 

From  fois,  it  is  clear  that  no  imaging  spectral  region  is  immune  from  contrast  reversal  effects.  It  is 
possible,  however,  to  reduce  their  magnitude,  or  in  some  cases  eliminate-  them  entirely  if  careful  nominal 
signature  prediction  is  utilised  together  with  criteria  for  eliminating  regions  where  large  signature  oscilla¬ 
tions  will  surely  exist.  A  more  detailed  discussion  of  this  problem  is  given  in  Section  3.  3. 

As  indicated  in  Figure  3,  a  reference  generation  proceas  is  used  to  develop  an  image  for  map-matching 
purposes.  Obviously,  to  ensure  systems  performance  fois  processing  stop  must  have  the  highest  dagraa  of 
accuracy  possible.  Two  types  of  predictive  •rrors  can  arive  from  lest  than  a  perfect  process.  The  first  ia 
the  result  of  having  to  synthetically  create  imagery  in  a  given  spectral  region  when  source  data  are  unavail¬ 
able.  The  second  involves  utilising  real  or  synthetic  reference  imagery  selected  or  generated  with  one  set 
of  envii'onmental  parameters  but  used  against  another  where  a  significant  signature  divergence  has  occurred 
The  mission  planner  should  use  a  nominal  rather  foan  abnormal  reference  image  when  large  signature  per¬ 
turbations  are  possible  which  can  not  be  accurately  predicted. 

When  spectral  band  conversion  is  necesaary  the  materials  within  foe  refere.ice  area  must  first  be  identi¬ 
fied.  The  synthetic  image  signature  is  generated  by  using  foe  known  physical  and  electrical  properties  of 
the  Identified  materials  in  conjunction  with  the  specified  Uluminator-target-detector  geometry. 

A  compilation  of  factors  influencing  foe  accuracy  of  reference  image  prediction  versua  the  actual  scene 
signature  is  presented  in  Table  6.  An  estimate  was  made  of  foe  expected  prediction  errora  for  homogeneous 
regiona  within  representative  reference  areas  for  each  spectral  region  and  is  given  In  Table  7.  Reaeonable 
uncertainty  values  of  perturbation  components  given  in  Tables  2,  4,  and  f>  were  used  tc  generate  foese  esti¬ 
mated  regional  errors.  Although  foese  values  should  only  be  used  as  a  guide,  they  can  provide  foe  mission 
planner  with  an  estimate  of  which  map-matching  algorifoms  can  not  ba  used  with  certain  forms  of  spectral 
imagery.  Thia  ia  due  to  foe  performance  breakdown  of  some  algorifom  classes  with  increased  regional 
errors.  The  estimated  regional  errors  in  Table  7  include  o>ntributions  from  material  identification  where 
appropriate. 

Results  given  in  Table  7  were  calculated  using  diurnal,  seasonal,  and  yearly  time-varying  signature  esti¬ 
mations  for  a  hypofoetical  reference  area  composed  of  45  paicent  vegetation,  30  percent  soil,  20  percent 
concrete,  and  S  percent  rock.  Snow/ice/water  complex  materials  were  excluded  from  fois  analysis. 
Vegetation  possesses  the  only  time-varying  dielectric  signature  variation  (excluding  snow/ice/water)  in  foe 
optical/uear  IR  region.  As  a  consequence  foe  error  bounds  given  in  Table  7  for  active  and  pasalve  types  in 
this  interval  should  be  evaluated  accordingly  when  ofoer  vegetation  proportiona  are  present.  Although  not  a 
factor  for  an  active  system,  large  actual  versus  predicted  error  bounds  for  passive  optical/near  IR  systema 
can  exist  if  diurnal  operation  is  desired  due  to  significant  spectral  irradlance  variations  present  in  day 
versus  ambient  night  light. 

As  in  foe  optical/near  IR  case,  foe  primary  aource  of  estimated  versus  actual  regional  error  bounds  in 
active  middle  and  thermal  IR,  and  microwave  images  ia  due  to  the  time-varying  vegetation  signature  pre¬ 
sent.  In  foese  Intervals,  however,  foe  general  lack  of  aource  data  necessitates  using  a  material 
identification  step  ia  producing  synthetic  reference  imagery.  The  resulting  errors  in  this  procedure  coupled 
with  the  lack  of  a  complete  physical  and  electrical  material  propertiea  catalog  produce  errors  in  foe  signa¬ 
ture  translation  process. 
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Tfcble  <>.  Parameter  Error  Impact  on  Intenaity  Eatimaie  Accuracy 
(Decreaalng  Order  of  Importance) 


IXES 

Paaaive 

Active 

Cptical/Near  IR 

Imaging  weathe  r 
alope/alope  aalmuth 

aeaaonal 

moiature  availability 
aurface  variationa 
dlu  rnal 

reflectance  knowledge 

Imaging  weathe  r 
slope /elope  aaimuth 

aeaaonal 

surface  variations 
moisture  availability 
reflectance  knowledge 

Middle  IR 

Imaging  weather 
thermal  inertia 
diurnal 

pre  -imaging  weather 
alope/alope  aaimuth 
moiature  availability 

aeaaonal 

aubaurface  variationa 
moiature  availability 
albedo  knowledge 
aurface  variationa 
emittance  knowledge 

Imaging  weather 
alope/alope  aaimuth 
aurface  variations 

seasonal 

reflective  knowledge 
moiature  availability 

Thermal  IR 

Imaging  weather 
thermal  inertia 
moiature  availability 
alope/alope  azimuth 
pre -imaging  weather 
aeaaonal 

aubaurface  variationa 
albedo  knowledge 
aurface  variationa 
diurnal 

emittance  knowledge 

Imaging  weather 
slope /slope  azimuth 
aurface  variations 
seasonal 

moisture  availability 
reflective  knowledge 

Microwave 

Moiature  availability 
thermal  inertia 
imaging  weather 
aurface  variationa 
emittance  knowledge 

Slope /slope  azimuth 
surface  variations 

seasonal 

moisture  availability 
reflectance  knowledge 

In  the  paaeive  middle  and  thermal  IR  caaea,  the  primary  aource  of  eatimated  veraua  actual  regional 
error  bounda  ia  due  to  the  ground  emiasion  component  governed  by  the  thermal  inertia  of  the  materiala  pre- 
aent.  In  addition  to  die  large  regional  error  preaent  between  moat  day /night  paira  analyaed  in  theae  caaea 
ia  die  fact  that  a  high  degree  of  anticorrelation,  indicative  of  the  inherent  contraat  reveraala,  alao  exlated. 
Theae  effecta  are  generally  noted  when  materi^a  with  a  moderate  to  high  range  of  thermal  inertiaa  are 
preaent  within  a  reference  area.  Aa  in  the  active  caaea  prevloualy  mentioned,  material  identification  errora 
and  data  gepa  in  material  property  librariea  alao  contribute  to  the  regional  errora  preaent. 

Aa  prevloualy  diacuaaed  in  thia  auction,  when  materiala  widi  high  microwave  emittancea  are  preaent 
widiin  a  reference  area,  the  reaulting  time -varying  paaaive  microwave  aignature  can  behave  aimilarly  to 
that  in  the  middle  and  thermal  IR  regie  '  .  If  materiala  with  noderate  to  low  microwave  emittancea  are  pre¬ 
aent,  the  variation  in  the  ground  temperature  component  of  the  apparent  brightneaa  temperature  due  to 
material  thermal  inertia  effecta  ia  damped,  and  a  greater  degree  of  regional  error  atabillty  reaulta. 

In  the  caae  of  range  or  phaa e-mod ulated  aenaora,  die  principal  aource  of  predicted  veraua  actual  regional 
errora  ia  due  to  the  time-varying  nature  of  vegetation  aignaturea  (particularly  deciduoua  treea)  preaent 
within  the  reference  area.  A  high  degree  of  reference  atabillty  ia  poaaible  widi  theae  aenaor  typea  if  care¬ 
ful  reference  aceno  acreenin;  ia  utiliaed. 
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Table  7  Estim/ted  Veraua  Actual  Regional  Error  Bound*. 


Expected  Error  Boards 

Spectral  Interval 

Type 

Low 

High 

Optical/Near  IR 

P 

16% 

25% 

A 

10% 

20% 

Middle  IR 

P 

20% 

>  100% 

A 

15% 

30% 

Thermal  IR 

P 

15% 

»  100% 

A 

10% 

30% 

Microwave 

P 

15% 

>100%*** 

A 

20% 

35% 

Range  or  Phase- 

A 

Small  except  when  deciduous  trees  present 

Mrdulated  Sensors 

*  Exclusive  cf  i'vc  w /ice /water  cjrnplex 


The  averag''.  tiij  foirenct  bet*>veen  the  ;.c.(.ual  mean  intensity 
level  difference  among  regions  and  the  predicted  intensity 
level  difference  among  regions  divided  by  the  actual 
intensity  range  period  spread  among  regions. 

*>!"«'  Whei.  vegetation  fraction  is  replaced  by  metal,  error 
bound  range  is  15%  to  30%. 


2,  ■‘i  Remedies 

Con*,  do.  reverssi'^.  nonsti uctured,  and  ^.rr.V-tive  errors  can  cause  map-m  tching  performance  degrade* 
tioni  '  an  other  i'>"'*or  types  (L  e. ,  i^eomet  distortions)  are  minimi.  There  are,  however,  four 

differ-  .it  r-  I'.cciy  cat.-,  xnat  can  minimize  the  impact  oi  these  errors  on  map-matching  systems 

performance.  They  i  dude:  (1)  accurate  nominal  signature  prediction,  (2)  proper  scene  selection, 

(J)  algorithm  flexibilsjy,  and  (4)  adaptive  performance  prediction. 

Accurate  nom-nal  signature  prediction  is  ^'.csirable  to  reduce  level  shifts,  hence  minimize  contrast 
reversals  between  homogeneous  regionc  -.^vithin  the  reference  area.  Errors  present  in  the  si(,nature  model, 
choice  of  nominal  atmospheric  conditions  or  material  identification  process  (if  utilized)  v>ill  all  contribute 
to  reduced  map-matching  performance.  Although  preprocessing  the  reference  and/or  sensor  scenes  can 
potentially  reduce  the  impact  of  global  and  local  bias  and  gain  changes,  the  results  are  quite  sensitive  to 
accurately  predicting  the  correct  time -varying  spatial  and  intensity  structure  of  the  image.  If  applied 
improperly,  preprocessing  steps  can  actually  reduce  rather  than  increase  systems  performance.  An  addi¬ 
tional  discussion  of  these  factors  is  given  in  Section  3. 

Proper  scene  selection  is  important  for  two  major  reasons.  Areas  that  are  prone  to  have  large  signature 
variations  in  a  given  spectral  region  due  to  contrast  reversals,  nonstructured  or  prediction  errors  should  e 
identified  and  eliminated  if  possible  in  the  scene  selection  process.  As  a  consequence,  an  accurate 
reference  scene  screening  procedure  is  desirable  so  that  an  estimate  of  the  probability  of  false  fix  (Pff)  can 
he  determined  under  a  variety  of  environmental  conditions  for  a  given  area.  It  is  necessary  here  to  evaluate 
-Ae  area  for  intrascene  redundancy  under  an  expected  operational  SNR.  If  an  unacceptably  high  Pff  results, 
'.'■'.e  cu.ididatc  reference  image  should  be  rejected.  A  more  detailed  pi.'esentation  of  these  topics  is  also 
given  in  Section  3, 

It  is  desirable  to  utilize  .ap-matching  algorithms  that  offer  a  degree  of  insensitivity  to  environmental 
perturbatloni:',  geometric  distortions  and  SNR  while. accurately  locating  the  true  match  point.  Each 
algorithm  cla  'correlation,  feature  extraction,  and  hybrid)  hax  its  own  advantages  and  disadvantages 
depending  on  the  type  of  inriagery  processed  and  the  magnitude  of  the  distortions  present.  A  more  detailed 
discussion  of  this  topic  is  given  in  Section  4. 

Since  map-matching  algorithm  performance  begins  to  break  down  with  increasing  distortion  present 
in  sensor  imagery,  it  is  desirable  to  utilize  a  technique  that  provides  a  confidence  estimate  of  the  quality  of 
the  fix.  A  generally  used  method  is  to  examine  the  surface  statistics  produced  by  the  map-matchi  -.g 
algorithm  correlation  of  reference  and  sensor  scenes.  Utilizing  a  simple  technique,  however,  that  does  not 
compensate  for  the  original  scene  properties  or  the  impact  of  the  algorithm  Itself  on  the  resulting  surface 
d.vtribution  can  be  inaccurate  when  typical  distortions  are  present.  A  more  detailed  discussion  of  adaptive 
performance  predictions  is  given  in  Section  5.2. 
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3.  Reference  map  construction  and  scene  selection  process 

Figure  4  describes  the  overall  map  construction  and  scene  selection  process.  Several  steps  are 
necessary  to  develop  a  reference  map  from  raw  seneur  data.  In  the  first  step  of  the  process,  it  may  be 
necessary  to  identify  the  scene  material  (especially  If  a  different  wavelength  Is  to  be  utilised)  and  geometri¬ 
cally  correct  for  other  viewing  aspects.  Once  this  is  accomplished  t'le  scene  signature  can  be  predicted  and 
a  nominal  signature  determined.  Due  to  environmental  factors  and  other  time-varying  variations  Inherent  In 
the  properties  of  scene.  It  is  also  necessary  to  predict  perturbations  to  the  nominal  that  are  likely  to  occur. 
Having  completed  the  signature  prediction  task  it  is  necessary  to  construct  th*  reference  map  and  check  via 
the  scene  selection  process  that  It  Is  adequate  for  the  map-matching  task.  In  constructing  the  reference  map 
in  many  cases  It  is  necessary  not  only  to  predict  intensity  levels  but  (depending  on  the  matching  algorithm) 
also  to  identify  homogeneous  regions  within  die  scene.  Once  this  is  accomplished  the  scene  can  be  checked 
via  mathematical  techniques  to  ensure  that  it  contains  sufficient  information  for  matching  purposes  and  that 
the  scene  is  sufficiently  unique  to  avoid  any  major  inter-scene  redundancy  problems.  Finally,  the  reference 
scene  must  be  tested  via  simulation  to  ensure  diat  it  is  suitable  under  re^  world  conditions. 

This  section  will  briefly  examine: 

1)  the  conversion  process, 

2)  the  problems  associated  with  signature  prediction, 

3)  construction  of  die  reference  map,  and 

4)  the  scene  selection  process. 


Conversion  Signature  Reference  Map  Scene 

Process  Prediction  Construction  Selection 


Figure  4.  Reference  map  construction  and  scene  selection  process 


3.  1  Conversion  process 


As  di  . cussed  previously  the  first  phase  of  reference  map  construction  generally  involves  conversion  of 
the  rav/  sensor  data:  1)  to  the  wavelength  or  frequency  of  the  sensor  onboard  the  vehicle  and  2)  to  the  geo¬ 
metrical  perspective  from  which  the  sensor  is  to  view  the  scene.  Because  the  raw  data  is  generally  not  at 
the  same  wavelength  as  the  sensor  it  is  necessary  to  estimate  the  material  properties  of  each  region  of  the 
scene.  Since  many  materials  have  very  similar  broad  band  reflecUnce  properties  in  the  optical/near  IR  por 
tion  of  the  spectrum  (from  which  most  raw  imagery  originates)  there  may  be  significant  mis -Identification 
errors  which  can  create  map-matching  difference  errors  and  ultimately  degrade  toul  system  performance. 

The  other  major  almost  InsurmounUbla  problem  is  to  adjust  the  imagery  for  the  geometry-  perspective 
which  the  sensor  is  likely  to  see.  For  systems  which  look  directly  down  (down-looking  systems)  the 
geometry  corrections  are  quite  simple  since  one  can  assume  a  flat  plane  model  for  the  ground.  For  other 
non-down  looking  systems  the  geometric  conversion  process  involves  developing  3-D  target  model  from  the 
original  2-D  imagery  and  then  creating  a  2-D  image  at  the  anticipated  perspective  angle.  Since  the  vehicle 
may  not  actually  fly  the  nominal  trajectory,  non  down-looking  systems  are  subject  to  geometric  errors 
which  require  significant  processing  efforts  to  remove. 


t nature  prediction 


Signatures  of  the  reference  map  need  to  be  determined  not  only  for  the  final  reference  map(s)  which  are 
stored  in  the  vehicle  for  comparison  but  also  to  test  (via  simulation)  the  performance  of  the  system. 
Seasonal  maps  of  t^'s  reference  area  may  need  to  be  developed  and  stored  for  use  in  at  least  some  map¬ 
matching  systems.  Depending  on  sensor  wavelength  and  map-matching  algorithm  it  may  also  be  necessary 
to  store  separate  ref>^rence  maps  which  account  for  diurnal,  atmospheric  and  meteorological  effects  on  the 
reference  scene  image.  The  mission  planner,  or  reference  scene  evaluator  is  cnutionsd  not  to  develop 
overly  sophistlcnted  tlgnnture  models  when  an  underspecified  set  of  conditions  will  exist.  Evan  worse  is 
the  case  where  poor  guesses  are  made  for  certain  input  variable  magnitudes;  since  in  some  cases  this  will 
result  in  nominal  reference  signature  with  greater  error  tiian  that  from  a  simple  model. 
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Perturbed  tignAture  v«rUtlont  from  the  nominal  ara  raqnirad  to  tost  the  performance  of  the  syetem 
under  a  variety  of  diurnal,  seasonal,  atmospheric,  and  meteorological  conditions.  One  should  utilise  this 
procedure  to  determine  whether  several  reference  maps  will  perform  better  under  a  variety  of  signature 
conditions  then  a  nominal  signature  prediction  which  is  not  accurate  for  any  one  scene  condition  but  is 
designated  for  compensating  for  these  variations!  effects. 

3. 3  Reference  man  construction 

In  the  reference  map  construction  area  there  a.v,  two  questions  which  need  to  be  addressed.  First,  what 
characteristics  should  the  reference  scene  possess?  Second,  how  should  the  reference  scene  be  evaluated? 

In  this  subsection  we  shall  briefly  discuss  the  choice  ol  a  reference  area.  In  the  subsequent  section  we  shall 
discuss  dve  simpler  question  of  reference  scene  evaluation. 

Table  8  lists  some  of  the  characteristics  in  the  ideal  reference  map  case  versus  the  real  world  situation, 
li  the  ideal  reference  map  characteristics  shown  in  this  table  existed  then  no  reference  screening  or  evalu¬ 
ation  procedure  would  be  necessary.  Philosophically,  one  can  not  control  mother  nature  nor  can  one  obtain 
agreement  (even  .f  one  could  control  mother  nature)  on  what  scene  characteristics  (number  of  homogeneous 
regions,  interpixel  correlation  length,  etc. )  are  best  for  map-matching  systems.  The  only  sure  thing  that 
can  be  said  about  reference  map  construction  is  that  certain  signature  characteristics  should  be  avoided,  and 
hence  this  is  the  major  topic  of  die  following  discussion.  Since  many  types  of  algorithms  require  that  homo¬ 
geneous  regions  or  features  be  identified  in  the  reference  map  a  brief  discussion  of  automatic  techniques  for 
the  region  extraction  is  included  here. 


Probable  Reference  Scene  Characteristics 
Probable  Case 


Table  8.  Ideal  Versus 

Ideal  Case 

Error  free  source  data  base 


No  reference  map  preparation  errors 


Reference  scene  contains 

-  A  single  homogeneous  region 

-  No  intra- scene  redundancy 

-  Statistically  independent  scene  elements 

-  Simple  statistical  intensity  distribution 

-  Time  and  space  invariant  signature 

without  contrast  reversals. 


Source  data  base  has: 

-  Finite  SNR 

-  Environmental  and  geometric  distortions 

present. 

Datum  plane  transferral  errors 
Imperfect  material  identification  and  signature 
models  used  in  spectral  translation. 

Reference  scene  usually  contains: 

-  Several  homogeneous  regions 

-  At  least  some  intra- scene  redundancy 

-  Interpoint  scene  element  correlation 

-  Complex  statistical  intenaity  distribution 
.  Time  and  space -varying  signature  with 

contrast  reversals. 


Proper  scene  selection.  Because  of  the  complexity  possessed  by  most  spectral  imagery  and  its  non¬ 
linear  space  and  time-varying  signature  characteristics,  the  reference  scene  selection  process  is  less  than 
an  exact  science.  It  is  generally  easier  to  make  qualitative  assessments  as  to  desirable  or  undesirable 
signature  physics  traits.  It  is  considerably  more  difficult,  however,  to  determine  exaefly  how  good  a  candi¬ 
date  reference  area  is  witiiout  rigorous  evaluation  due  to  (he  statistical  nature  of  expected  environmental  and 
geometric  distortions,  SNR  effects  and  intrascene  redundancy. 

The  net  effect  of  these  degradations  is  to  impact  map-matching  algorithm  performance,  and  hence,  die 
reliability  of  the  fixing  process  Itself.  An  examination  of  algorithm  class  sensitivity  to  SNR  and  contrast 
reversals  is  given  in  Section  4.  3,  while  a  review  of  fix  performance  estimation  measures  is  presented  in 
Section  5. 2. 

If  a  map-matching  algorithm  is  utilised  which  ie  sensitive  to  contrast  reversals  (i.  e. ,  ordinary  correla¬ 
tion  metrics)  then  vegetation  diat  exhibits  strong  time-varying  growth  variations  should  be  omitted  or 
minimised  in  update  areas  in  every  spectral  interval.  Similarly,  it  is  also  advisable  tc  eliminate  candidate 
update  areas  where  low  material 'thermal  inertia  and  short  wavelength  reflectance  in  the  passive  middle  and 
thermal  IR,  and  microwave  (for  high  microwave  omittance  materials)  regions  predominates  to  avoid  contrast 
reversal  effects.  From  Table  5,  it  is  clear  that  the  snow /ice /water  complex  can  adversely  alter  the 
reference  area  signature  in  each  epeotral  interval.  Obviously  then,  water  bodies  should  only  be  included  in 
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reference  area*  if  they  are  unlikely  to  freeae  becauae  of  the  moderate  to  large  algnature  pertnrbr.tioni  that 
can  reault  in  active  and  paaeive  apectral  imagery.  Unleaa  phaae -modulated  or  range>gated  ayatema  are 
utilik:ed,  diaaatroua  fixing  reaulta  will  often  occur  with  ordinary  correlation  or  feature  matchinf  algoridima 
when  anow /ice/water  ia  preaent  and  changea  in  complex  atate  art  expected. 

If  map-matching  algo  rithma  are  uaed  which  are  aenaitlve  to  SNR  (primarUy  feature  matching  and  to  a 
leaaer  degree  die  hybrid  proceaaing  approach),  then  ragiona^ere  atrong  at^-:apheric  and  meteorological 
variationa  exist  ahould  be  carefully  evaluated.  The  impact  of  atmoapheric  parametera  (particularly 
attenuation  and  aeroaol  acattering)  typically  decreaaea  with  increaaing  wavelength,  but  atlll  generally  forma 
the  limiting  operational  caae  to  the  thermal  IR  region  (where  reradiation  becomea  important).  In  the 
passive  microwave  region,  reradiation  from  percipitabla  water  can  introduce  small  to  large  signature  vari¬ 
ations;  partioilarly  when  materials  of  differing  microwcve  amittancea  exist.  Radars,  however,  are 
generally  unaffected  by  all  but  the  modt  severe  atmot  pher'.c  conditions. 

Although  meteorological  effects  (i.  e. ,  latent  and  sensible  heat  transfer)  typically  produce  a  entailer  per¬ 
formance  degradation  than  atmoapheric  ones,  they  directly  impact  the  terrain  signature  in  each  spectral 
region  when  soil  moisture  ia  present  by  governing  its  rate  of  avaprration.  For  each  active  spectral  region 
and  passive  optical/near  IR,  this  appears  aa  a  change  in  soil  reflectance.  In  the  passive  middle  and  thermal 
IR,  and  microwave  regioas  soil  moisture  variations  alter  the  emissive  powers  of  the  surface. 

Soil  moisture  effects  will  gunerally  impact  or  ;;inary  correlation  and  feature  matching  algorithm  per¬ 
formance  the  greatest,  since  its  presence  rn  sensor  imagery  is  space  and  time-varying  and  is  often  not 
equally  proportioned  widiin  a  homogeneous  region.  The  impact  of  latent  and  sensible  heat  transfer  for  low 
soil  moisture  and  high  atmospheric  precipitable  water  will  generally  be  to  reduce  the  dynamic  range,  hence 
contrast  between  homogeneous  regions,  in  passive  middle  and  thermal  IR,  and  microwave  imagery. 

In  some  cases  even  these  prccedures  'vill  be  inadequate  to  produce  representative  Imagery  for  guidance 
updating  purpose*.  Here,  it  may  be  necessary  to  select  multiple  reference  images  of  the  same  area  to 
compensate  for  diurnal  and  seasonal  effects.  From  this,  ^e  mission  planner  can  either  select  the  niost 
representative  image  in  near  real-time  or  store  the  set  of  multiple  frames  in  the  vehicle. 

Diurnal  variations  in  passive  middle  and  diermal  IR,  and  microwave  imagery  tend  to  be  region-based. 
Seasonal  variationa  except  those  induced  by  snow /ice /water  tend  to  be  interregion.tl  for  all  the  candidate 
sensor  types  considered  here.  As  a  consequence,  the  hybrid  nnap-nuitching  algorldtm  is  often  desirable  if 
an  adequate  SNR  exists.  From  this,  it  is  apparent  that  proper  algorithm  selectior  for  a  given  sensor  type 
can  often  simplify  the  task  of  nominal  reference  scene  prediction.  Conversely,  r.siog  a  sub-optimal 
algorithm  will  often  place  an  overly  stringent  accuracy  requirement  on  signature  prediction,  and  signifi¬ 
cantly  increase  the  time  required  for  reference  scene  preparation. 

Preprocessing  references  and  sensor  images  or  using  binary  data  correlation  can  reduce  the  impact  of 
signature  perturbation  factors.  As  previously  discussed,  such  schemes  can  n<.t  bo  successfully  utilized 
without  a  borough  understanding  of  the  expected  imaging  physics,  SNR  and  geometric  distortions  present. 

If  applied  blindly,  these  techniques  can  often  reduce,  rafter  than  enhance,  guidance  updating  systems  per¬ 
formance. 

Because  of  the  inherent  deficiencies  in  nominal  signature  prediction  for  a  given  rensor  type  coupled  wift 
map-matching  algorithm  limitations,  it  is  often  desirable  to  employ  adaptive  performance  prediction 
measures  to  estimate  the  quality  of  Individual  fixes,  A  discussion  of  possible  performance  prediction 
techniques  for  guidance  updating  applications  and  ften  limitations  is  given  in  Section  5.  Z. 

Region  extraction  (12-22).  Obviously  homogeneous  regions  or  features  in  fte  scene  cau  be  found  visually; 
however,  when  scenes  are  described  digitally  by  large  arrays  of  numbers,  it  is  highly  desirable  to  intro¬ 
duce  some  level  of  automation  into  the  process.  There  are  two  different  approaches  for  automatically 
extracting  regions  from  scenes:  1)  those  based  on  edge  operators  and  2)  ftose  based  on  the  stationarity  pro¬ 
perties  of  the  region. 

Fdge  approaches  apply  gradient  or  Laplacian-type  operators  to  the  scene  and  then  use  threshold  techni¬ 
ques  to  decide  upon  fte  existence  of  any  edge  (the  boundary  of  a  feature).  The  major  danger  in  using  these 
techniques  is  that  noise  and  distortion  can  make  it  very  difficult  to  locate  edges  in  the  sensor  innagery. 

Homogeneous  regions  may  also  be  located  using  the  statistical  property  of  stationarity  (first  order,  con¬ 
stant  mean  level  in  the  region;  second-order,  mean  and  variance  constant  and  autocorrelation  independent  of 
position).  In  this  area-based  approach,  one  would  grow  regions  of  spatially  connected  pixels  which  would 
possess  ftis  property.  While  this  approach  is  less  susceptible  to  problems  of  noise  and  distortions  it  is 
computationally  more  complex  than  the  simpler  edge  operator  techniques. 

3.4  Scene  selection 

The  scene  selection  process  is  concerned  with  choosing  maps  for  which  the  probability  of  matching  a 
sensor  image  from  within  the  reference  map  boundary  is  high.  This  process  has  both  physical  and 
mathematical  implications.  There  will  obviously  be  signature  differences  between  the  sensed  innags  and  its 
reference  map  counterpart  due  to  such  factors  as  geometric,  atmospheric,  meteorological,  diurnal,  and 


•eascnal  effects.  These  effect*  on  tystein  perform&nc»  can  be  tnu.^'nlsed  In  the  extreme  by  either  pre¬ 
paring  the  ref<irence  map  to  be  near  real-time  estimate  of  the  sensor  image  at  ttie  moment  iv  overflies  the 
reference  area  or  by  developing  scene  particular  algorithms  that  are  relatively  invariant  to  the  signature 
deviations  between  the  sensor  image  and  reference  map.  Realistically  one  must  reach  a  comp.-omise 
between  these  two  extremes  and  develop  a  reference  map  which  will  reduce  the  signature  deviations 
especially  in  defining  the  boundaries  of  a  homogeneous  region  and  utilise  a  matching  algoridim  that  will 
compensate  for  any  remaining  signature  differences  between  the  two  maps. 

In  general,  successive  screening  techniques  from  eimple  math  tests  to  full-blown  simulation*  are 
chosen  and  used  to  evaluate  the  candidt  ce  reference  area.  Since  computer  processing  requirements 
increase  consid<*ri.'jly  with  each  screening  *i'.ep,  it  is  desirable  f(>r  unacceptable  reference  scenes  to  be 
identified  and  rejected  before  the  final  simulation  analysis  if  possible. 

The  mathematical  criteria  for  reference  scene  selection  requires  that  there  be  (1)  sufficient  information 
for  map'-inetching  and  (2)  a  minimal  amount  of  Intrascene  spatial  redundancy  within  ^e  reference  map 
boundary.  Techniques  exist  for  measuring  the  information  content  in  the  scene  to  ensure  that  the  sensor 
image  size  (In  terms  of  resolution  elements)  contain  a  sufficient  number  of  independent  elements.  The  more 
difficult  issue  and  yet  unresolved  is  the  determination  of  measure  for  scene  uniqueness.  Equipped  with  such 
a  measure  it  would  be  poscible  within  a  reference  map  boundary  to  test  the  ensemble  of  possible  sensor 
images  to  determine  the  amount  of  intrascene  spatial  redundancy. 

Heights  of  the  secondary  correlation  peaks  and  their  nr\agnitude  determined  by  autocorrelating  a  particu¬ 
lar  sensor  map  over  the  reference  map  area  yield  the  location  of  areas  where  there  is  a  possible  spatial 
redundancy  problem.  Two  problems  emerge  from  attempting  to  use  this  as  a  measure  of  the  uniqueness 
problem.  First,  in  real  world  imagery  the  magnitude  of  &e  intensity  level  shifts  witiiir  the  imagery  may 
have  a  significant  impact  on  the  location  of  secondary  peaks.  Thus  this  approach  does  not  seem  fruitful  for 
measuring  scene  uniqueness.  Second,  this  approach  use*  texture  information  within  a  region  which  may  or 
may  not  be  used  in  the  matching  algorithm;  consequently,  the  results  may  be  different  when  texture  informa¬ 
tion  is  omitted. 

The  underlying  spatial  patterns  in  the  map  as  designated  by  the  sise  and  sh  >pe  of  the  homogeneous 
regions  are  the  primary  concern  in  dealing  with  the  spatial  redundancy  probl«..'n.  One  method  for  measuring 
scene  uniqueness  would  be  to  use  the  correlation  surface  associated  with  a  specialized  hybrid  algorithm  4* 
a  mean*  for  screening  reference  maps.  Here  the  reference  area  would  be  broken  up  into  homogeneous 
regie  \*  and  each  pixel  within  the  region  would  be  replaced  by  the  mean  intent 'ty  level  of  the  entire  region. 

An  autocorrelation  of  a  particular  sensor  map  over  the  reference  map  would  be  performed  using  a  hybrid 
correlation  algorithm.  High  secondary  correlation  peaks  would  indicate  area*  where  spatial  scene 
redundancy  potentially  could  be  a  problem.  By  pulling  out  a  number  of  sensor  maps  from  the  reference  map 
boundary  and  repeating  this  process  one  could  determine  as  first-order  measure  ^e  uniqueness  properties 
of  the  reference  map. 

The  most  accurate  evaluation  procedure  uses  a  Monte  Carlo  simulation  lo  provide  an  estimate  of  the 
update  circular  error  probability  (CEP)  and  Pff  for  a  given  reference  area  under  a  variety  of  condition*  (23). 
Samples  are  drawn  from  statistic^  distributions  that  represent  vehicle  position-velocity  characteristics, 
and  are  used  to  specify  the  sensor  scene  location  within  the  reference  image.  Samples  from  another  distri¬ 
bution  are  used  to  specify  the  imaging  environmental  characteristics  present  (i.  e.,  time  of  day).  An 
intensity  computation  for  the  specified  conditions  is  performed  for  each  subscene  location  by  using  the 
appropriate  signature  model.  Noise  terms  and  geometric  distortions  are  similarly  introduced  into  the 
sensor  scene. 

The  map-matching  operation  is  dien  performed  between  the  reference  and  specified  synthetic  sensor  maps 
using  the  selected  matching  algorithm.  The  along  and  cross-track  difference  between  the  computed  and 
actual  (sensor  scent'  draw)  match  point  locations  is  determined,  and  by  using  a  predetermined  criteria,  the 
update  is  catalogued  a*  eitlier  a  good  or  false  fix. 

Since  each  of  the  variables  are  represented  by  statletical  distributions,  the  simulation  can  be  run  a 
specified  number  of  times  to  provide  CEP  and  Pff  estimates  over  the  range  of  expected  update  conditions. 
From  this,  the  reference  scene  suitability  for  guidance  updating  applications  can  be  determined  versus  a 
predetermined  performance  criteria. 

An  estimate  of  the  spatial  redundancy  present  within  the  reference  scene  is  provided  by  this  procedure 
since  the  location  of  the  sensor  scene  is  randomly  selected  from  within  the  reference  map  boundaries.  Like¬ 
wise,  estimate*  of  the  Impact  of  environmental.  SNR  effect*  on  reference  mr.p  performance 

are  also  evaluated  by  this  procedure. 

It  should  be  recognized,  however,  that  the  quality  of  die  performance  estimate  provided  by  a  Monte  Carlo 
simulation  for  a  given  reference  area  is  generally  a  strong  function  of  die  preprocessing  and  map-matching 
algoridims  and  the  characteristics  of  the  environmental,  and  gsometric  distortions  and  SNR  selected.  The 
use  of  different  preprocessing  or  map-matching  algovithms  for  refersnee  scene  evaluation  and  guidance 
upd'^.ting  should  be  avoided  to  minimise  performancs  degradations. 
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Any  uncerUinty  In  ■pecifyins  the  random  variable  dietribution  propertiea  utUlced  in  the  eimolation  will 
reiult  in  a  decreased  confidence  in  the  reference  ecene  evaluation  proceaa.  If  large  uncertainties  exist  in  a 
given  variable  distribution,  it  is  better  to  eliminate  the  variable  from  the  simulation.  If  significant 
uncertainty  errors  are  present  in  several  distributions,  the  benefits  of  employing  this  form  of  reference 
Scene  evaluation  decrease  and  the  resulting  performance  estimates  produced  are  often  unreliable. 


4.  Alttoriihm  selection  process 

Figure  S  describes  the  overall  algorithm  selection  process.  The  sensor  image  io  generally  considered 
to  be  some  subset  of  the  reference  map  corrupted  by  errors.  A  matching  algorithm  is  then  used  to  locate 
the  position  Jt  the  sensed  image  in  the  reference  map  coordinate  system.  Based  on  an  analysis  of  system 
performance,  it  is  possible  to  determine  the  capability  of  each  algorithm  to  accommodate  various  types  of 

TOT.  Ultimately,  since  for  each  sensor  type  some  errors  are  more  dominant  than  others,  it  is  possible 
to  determine  the  most  appropriate  algorithm  for  each  sensor  type. 

This  section  will  discuss  the  following  topics:  (1)  a  description  and  categorisation  of  error  sources; 

(2)  u  description  and  classification  of  matching  algorithms;  and  (3)  an  analysis  of  the  compatibility  of  various 
Ugorlthms  to  accommodate  different  error  sources. 

4,  1  Error  sources 

The  problems  associated  with  image  dynamics,  geometrical  distortions,  noise,  and  other  error  sources 
can  be  lumped  into  four  mutually  exclusive  comprehensive  categories.  These  categories  are  defined  as; 

1)  Global  errors--those  errors  which  uniformly  affect  the  intensity  level  of  all  scene  elements.  This 
category  would  include  geometric  distortions  and  bias  and  gain  changes. 

2)  Regional  errors--those  errors  where  the  charts  in  intensity  levels  occurs  uniformly  only  within 
homogeneous  regions  within  the  scene.  Examples  would  include  region*level  shifts  (contrast 
reversals)  due  to  image  dynamics  and  predictive  coding  errors  from  incorrect  reference  map 
construction. 

3)  Local  errors  "-errors  expected  to  affect  each  pixel  or  grouping  of  pixels  (contained  within  an 
interpixel  correlation  length)  independently.  The  primary  example  of  this  error  source  is 
additive  noise. 

4)  Nonstructured  errors- -this  is  a  catchall  category  designed  to  fit  those  errors  whose  effect  on  the 
scene  can  not  be  described  as  being  global,  regional,  or  local  (an  example  being  a  cloud  cover  over 
portions  of  the  target  area  changing  the  signature  in  a  nonstructured  manner). 

The  advantage  of  this  formulation  of  the  error  source  is  that  by  grouping  errors  into  these  categories 
it  is  simpler  to  describe  the  utility  of  each  algorithm  relative  to  a  given  class  of  error  source  radier  than 
having  to  backtrack  and  deal  with  each  error/algorithm  combimaion  on  an  individual  basis. 


Reference  Image 


Appropriate 
Algo  rlthm/ 
Sensor 
Combinations 


4.  2  Map-Matching  Aliorlthm  Classes 


There  are  three  classes  of  algorithms  which  can  be  employed  to  perform  the  image  matching  task.  These 
algorithms  include  correlation,  feature  matching,  and  hybrid  classes  (24). 
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All  algorithms  perform  four  operations;  (1)  transfornnation  of  the  original  Intensity  data  associated  with 
each  resolution  element  In  u&cn  sensed  Image  and  reference  map;  (2)  establishment  of  a  metric  for  com¬ 
paring  a  portion  of  the  reference  map  to  the  sensed  image;  (3)  the  computation  of  that  metric  for  all  possible 
positions  of  comparison  between  the  reference  and  sensor  maps;  and  (-1)  a  selection  rule  (generally  the 
extrennum  metric  value)  for  delineating  the  match  position  baaed  on  the  metric  value. 

Correlation  types  of  algorithms  use  the  intensity  values  associated  with  the  resolution  elements  of  each 
map  (of  aome  transformation  of  these  intensity  values,  i.e. ,  normaliaation)  as  the  map  data  to  be  used  in 
computing  the  metric.  Correlation  metrics  measure  either  the  degree  of  similaritv  (i.  e. ,  product  type 
algorithm)  or  dissimilarity  (i.e. .  difference  squared)  between  the  sensed  image  and  the  portion  of  the 
reference  map  it  it  being  compared  against. 

Feature  matching  algorithms  do  not  utilise  intensity  data  per  se  but  attempt  to  work  with  only  features  in 
the  scene  (25).  This  is  generally  accomplished  by  using  algorithms  to  locate  boundaries  or  edges  between 
regions.  Edge  or  boundary  Information  is  extended  to  determine  the  position  at  which  boundaries  or  edges 
intersect.  The  position  of  this  vertex  point  and  the  direction  and  number  of  line  segments  emanating  from 
the  vertex  point  form  the  basis  for  map  comparisons  wldi  the  metric  being  some  form  of  a  mean  square 
distance  measure  between  locations  of  vertices  in  the  reference  and  sensed  map.  This  distance  measure 
may  be  weighted  by  the  number  and  direction  of  line  segments  emanating  from  the  vertex  point  with  multiple 
intersecting  vertices  being  weighted  more  heavily. 

The  hybrid  algorithms  (26)  uses  a  combination  of  Intensity  level  and  region  identificatir n  information  in 
determining  a  match  location.  In  this  class  of  algorithm  homogeneous  regions  in  the  ref'jrence  scene  are 
identified  and  all  resolution  elements  within  the  region  are  tagged  as  belonging  to  the  region.  When  the 
sensor  is  compared  to  a  portion  of  the  reference  map,  an  assumption  is  made  that  this  position  of  compari¬ 
son  is  the  correct  one,  and  the  sensor  image  is  broken  up  into  homogeneous  regions  as  identified  Sy  the 
counterpart  elements  of  the  reference  map.  The  elements  in  each  region  are  correlated  separately  using  a 
correlation  algorithm,  and  the  total  correlation  between  the  two  maps  is  found  by  summing  the  individual 
correlations  taken  over  each  homogeneous  segment  of  the  reference  map. 

4.  3  Analysis  of  algorithm  compatibility 

Let  us  consider  which  class  of  algorithm  is  most  appropriate  for  accommodating  each  class  of  error 
source  separately.  Table  9  shows  a  rating  of  the  algorithm's  ability  to  accommodate  each  error  class. 
Examining  the  errors  relative  to  the  algorithms,  all  algorithms  can  readily  accommodate  global  errors. 
Correlation  and  hybrid  algorithms,  however,  probably  have  somewhat  more  difficulty  in  accommodating 
this  type  of  error.  Corrective  action  for  compensating  for  global  errors  nclude  processing  of  sensor  ele- 
ments  in  smaller  spatial  groups  to  accommodate  geometric  errors  (27-29),  normalisation  of  intensity  levels 
to  compensate  for  bias  errors  and  gain  shifts,  and  extending  the  algorithm  ,  earch  dimension  to  include 
searching  the  scene  for  scale  and  rotational  effects.  Correlation  and  hybrid  algorithm  corrective  measures 
would  rely  most  heavily  on  spatial  grouping  and  intensicy  level  compensation,  while  feature  matching  algo¬ 
rithms  (working  with  It.'ss  data  to  begin  with)  would  primarily  resort  to  search  techniques  to  compensate  for 
global  errors. 


Table  9.  Algorithm  Ability  to  AccoRunodate  Each  Class  of  Error 

Error  Class 


Algorithm 

Global 

Regional 

Local 

Nonstructured 

Correlation 

Medium 

Poor 

Good 

Good 

Feature 

Matching 

Good 

Good 

Poor 

Poor 

Hybrid 

Medium 

Good 

Medium 

Good 

Correlation  algorithms  are  extremely  poor  performers  in  the  pn  sence  of  regional  errors,  the  possible 
solution  being  (besides  switching  to  one  of  the  other  algorithms)  to  stoye  and  search  over  multiple  reference 
maps,  restrict  the  wavelength  of  the  imagery  to  spectral  regions  where  regional  errors  are  not  likely  to 
occur,  or  to  locate  reference  maps  in  geographic^  areas  in  which  regional  errors  are  unlikely.  Both  the 
feature  matching  and  hybrid  class  of  algorithms  are  good  in  accommodating  regional  errors. 

Local  errors  such  as  noise  can  cause  significant  problems  in  the  performance  of  feature  matching 
algorithms  primarily  due  to  the  difficulty  in  extracting  features  of  lino  boundaries  from  dte  sensed  imagery 
using  edge  type  operators.  Correlation  type  algorithms  are  virtually  immune  to  local  errors,  while  hybrid 
algorithms  are  susceptible  to  this  error  source  if  there  is  also  a  scene  redundancy  problem  with  noise, 
making  it  more  difficult  to  distinguish  images  with  similar  spatial  patterns.  The  only  corrective  meaaure 
for  feature  matching  algorithms  in  the  presence  of  local  errors  is  to  switch  to  one  of  die  other  two  classes 
of  algorithms. 


Finally,  since  feature  matching  algorithms  use  less  information  in  tiie  scene  than  the  other  two  types  of 
algo  ri^hme,  they  are  most  susceptible  to  nonet ructured  errurs  where  positions  of  the  sensed  image  may  look 
oblHerated  when  compared  to  the  reference  map.  Correlation  and  hybrid  algorithm. i  can  still  perform  quite 
well  even  in  the  presence  of  large  miasi->g  areas  in  the  sensed  image. 

A»  discussed  above,  each  algorithm  has  advantages  and  disadvantages  relative  to  certain  types  of  erroi 
sources;  however,  real  world  systems  are  likely  to  be  faced  with  a  combir.ii  ?n  of  error  sources  to  deal 
with. 

From  the  discussion  In  Section  2,  it  Is  seen  that  certain  sensor  bands  have  characteristic  errors  pri¬ 
marily  regional  (i.  e. ,  contrast  reversals  and  predictive)  errors  and  local  ei  rors.  i***  magnitude 

of  these  '.rrors  sources  (and  excluding  the  effects  of  global  and  nonstructuru.i  errors;  it  is  possible  to  deter¬ 
mine  the  compatibility  of  sensors  with  matching  algorithms.  If  regional  ei.-crs  donninate  the  process,  then 
a  feature  matching  algorithm  is  most  appropriate.  If  local  errors  dominate,  t.?en  correlation  algorithms 
look  most  attractive.  In  the  presence  of  both  local  and  regional  errors  then  the  >  rid  class  of  algorithm 
looks  most  appropriate. 

To  summarise,  all  error  sources  can  be  placed  into  one  of  four  generic  categories.  Using  these  cate¬ 
gories  one  can  analyze  the  performance  of  dte  three  types  of  algorithms  relative  to  a  particular  error 
source.  Some  algorithms  are  more  capable  than  others  at  accommodating  certain  classes  of  error.  In  the 
end  the  final  algorithm  choice  wUl  depend  os.  determining  the  weighting  of  the  error  sources  that  the  system 
is  likely  to  encounter. 

An  analysis  was  performed  to  determine  die  optimum  map-matching  algorithm  class  for  each  sensor 
operating  band  based  on  the  regional  errors  given  in  Table  7,  as  well  as  sensor  characteristics  and  oper¬ 
ational  considerations.  The  results  of  this  analysis  are  summarized  in  Table  10.  Aldiough  in  no  sepsor 
case  is  one  algorithm  class  superior  to  the  others,  a  number  of  caveats  have  been  developed  and  presented 
as  a  guide  to  the  mission  planner  to  ensure  optimum  performance. 


Table  10.  Map>matching  algorithm  selection  based  on  designated  sensor  operating  region 


Sensor  Region 


lYpe 


Algorithm  Selection* 


Optical /Near  IR  P  It  A 


Middle  m 


P 


A 


•  Correlation  when  SNR  low 

-  Hybrid  when  SNR  moderate  and  vegetation  present 

-  Feature  matching  w’hen  SNR  high  and  vegetation  absent 

-  Correlaticn  unacceptable  because  of  regional  thermal 
inertia  effects. 

-  Hybrid  when  SNR  low  to  moderate  and  vegetation  pres'^nt 

-  Feature  matching  generally  undesirable  unless  high  SNR 
exists  and  vegetation  is  absent. 

•  Same  as  Optical /Near  IR 


Thermal  IR 


P 


-  Same  as  Passive  Middle  IR 


Microwave 


A 

P 


A 


Range  or  Phase  -  A 

Modulated  Sensors 


-  Same  as  Optical/Near  IR 

-  Correlation  when  SNR  low  and  microwave  reflectance 
dominates. 

>  Hybrid  when  SNR  moderate  and  vegetation  present 

-  Feature  matching  when  SNR  high  and  vegetation  absent 

-  Correlation  when  SNR  low  or  with  one -dimensional  imaging 

systems. 

-  Hybrid  when  SNR  moderate,  vegetation  present  or  with  two- 
dimensional  imaging  systems. 

-  Feature  matching  generally  unacceptable  because  of  inade¬ 
quate  SNR  unless  specialized  preprocessing  used. 

-  Correlation  only  when  Tow  SNF  present 

-  Hybrid  unnecessary  since  regional  er^'ors  are  generally 
small. 

-  Feature  matching  desirable  when  SI'S!  high 


*  When  moderate  nonatructured  errors  or  snow/ice/water  are  present,  tite  hybrid  approach  must  be  'taed 
with  all  systems  except  range  or  phase-modulated  sensers  to  ensure  update  reliability. 
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CorrelaHon  It  detirtble  in  activ*  and  pattlve  optical/near  IR  catet  where  a  low  SNR  it  pretent.  Sourcec 
of  thlt  type  Include  pattlve  low-light  level  operation,  low  acene  contraat  operation,  or  when  a  high  atmoe- 
pheric  aeroeol  content  It  present  In  the  Imaging  path.  When  eigniflcant  vegeUtlon  la  preeent  in  die 
reference  area,  the  hybrid  approach  la  detlrable,  whUe  feature  extraction  ie  reaerved  for  caeea  when  a  high 
SNR  nxlatt  and  vegetation  diat  exhlbita  a  time-varying  regional  boundary  ahift  la  abaert. 

Correlation  claaa  algorlthma  provide  unacceptable  performance  with  paaelve  middle  and  thermal  IR  and 
generally  with  paaaive  microwave  imagery  (when  microwave  emlttance  predoiriiiwiteB)  due  to  contraat 
revertalt  or  moderate  to  large  time-varying  regional  errora  induced  by  materia',  thermal  Inertia  effectt. 
Traditional  correlation  technlquea, Including  binary  converalon  preproceaalng,  gen*  rally  provide  Inalgnifl- 
cant  performance  improvement  when  applied  to  middle  and  thermal  IH  Imagery  and  amall  to  moderate 
improvem<'nta  with  paaelve  microwave  Imagery.  Th«  hybrid  approach  in  preferable  in  caaea  where  lea* 
than  a  high  SNR  exinta  due  to  aenaor,  imaging  contraat  or  atmoapheric  reradlatlon  conaideratlona,  or  when 
a  time-varying  vegetation  algnature  la  preeent.  Feature  matchi^  application  la  operationally  limited  aa  In 
other  apectral  regiona  to  caaea  where  a  high  SNR  exiata  and  time- varying  vegetation  coverage  ia  abaent. 

Commenta  given  for  the  optical /near  IR  region  are  generally  applicable  for  all  map-matching  ayatema 
uaing  active  aeneora.  The  principal  limitation*  of  active  middle  and  (particularly)  titermal  IR  ayatema  for 
epplicatlona  againat  natural  materiala  1*  ttie  low  imaging  contraat  typically  preaent.  It  la  often  neceaaary  to 
utUize  dynamic  range  expanaion  preproceaalng  technlquea  with  theae  aenaor  type*,  which  limit*  the  uae  of 
feature  matching  metiiod*  in  theae  caaea  unleaa  a  high  data  SNR  exiata. 

Aldtough  atmoapheric  effect*  generally  have  a  negligible  Impact  on  radar  image  contraat,  the  moderate 
to  low  SNR  typically  preaent  for  moat  propoaed  miaai!  j-bome  ayatema  coupled  wlti»  apecular  material 
return*  fro.T>  the  reference  area  provide  other  problem*  for  operational  map-matching  ayatema.  With  one- 
dlmenalonal  radar  map-match:  ig  ayatema  correlation  algorithm*  are  uaually  preferred  over  feature 
matching  to  minimlxe  the  number  of  dlacrets  acattera  required  to  enaure  adequate  performance.  Hybrid 
algorithm*  are  preferable  when  a  moderate  SNR  exiata  or  when  aignlllcant  vegetation  la  preaent  that 
poaaeaaea  a  time-varying  algnature.  Feature  matching  algorithm*  are  generally  unacceptable  lor  pro- 
ceaaing  miaalle  borne  radar  data  becauae  of  typically  low  SNR*  unleaa  apeclallaed  preproceaalng  technlquea 
•re  u*ed  which  emphaaiae  atable, while  deemphaaialng  potentially  un*table,edgea  preaent. 

For  range  or  phaae-modulated  aenaor*,  the  hybrid  approach  1*  generally  unwarranted  (unleaa  a  aignlfl- 
cant  amount  of  deciduou*  tree*  exiat)  becauae  of  the  generally  time -Invariant  nature  of  tiieae  form*  of 
reference  imagery.  The  choice  between  no.  relation  and  feature  matching  approache*  here  ahould  be  deter¬ 
mined  veraua  the  expected  SNR  aince  predictive  and  nonatructured  errora  are  generally  email. 

In  addition  to  the  caveata  juat  presented.  It  ahotild  be  recognlaed  that  other  error  type*  aometlme*  pre¬ 
eent  can  aigniflcantly  alter  map-matching  performance.  Correlation  and  feature  matching  algorltiim  claaa 
performance  are  aenaltlve  to  predictive  (U  e. ,  anow/lce/water  complex)  and  nonatructured  (1.  e. ,  ahadowing) 
reference  map  error*.  In  the  event  that  a  high  probabUlty  of  time  and  apace-varying  anow/lce/water  or 
ahadowing  exiat*  within  the  reference  area,  the  hybrid  algorithm  claaa  1*  preferable.  The  only  practical 
exception  to  thla,  ellowing  adequate  correlation  or  feature  matching  performance,  would  Involve  die  blockage 
of  only  a  email  amount  of  the  total  map  Infornnatlon  content  (i.  e. ,  number  of  independent  element*)  tnd/or 
total  map  area. 

In  aummary,  when  aenaor  character!* tica  or  operational  con* ide ration*  reault  In  a  low  SNR  and  when  die 
aelected  reference  area  ha*  a  relatively  time  invariant  algnature,  correlation  claaa  algorithm*  ahould  be 
conaidered.  When  a  moderate  SNR  exiata  and  predictive  and  nonatructured  errora  are  expected,  the  hybrid 
claaa  ia  preferable.  In  caaea  where  a  high  SNR  exiat*  and  predictive  and  nonatructured  errora  are  email, 
feature  mapping  1*  deairable. 


Performance  prediction 


Major  performance  conaideratlona  for  image  matching  ayatema  involve  (1)  the  avoidance  of  falae  flxea 
during  acquieition,  and  (2)  the  accuracy  with  which  die  poaition  fix  can  be  mad*.  The  major  focua  of  thi* 
paper  i*  on  the  acquiaition  phaa*  of  image  matching,  which  i*  the  more  difficult  and  important  part  of  the 
overall  problem  to  be  aolved.  The  acquiaition  ayatem  deeigner  relative  to  performance  meaaurea  ia  con¬ 
cerned  (1)  with  developing  general  guideline*  for  performance  aa  a  function  of  aenaor  and  computational 
algorithm  characteriatlca,  and  (2)  real-time  acene  dependent  eatimatea  of  ayatem  performance  in  order  to 
determine  whether  or  not  a  poaition  fix  ia  valid.  Thi*  aection  ia  therefore  divided  into  two  parte;  one 
dealing  with  the  general  development  of  performance  guideline*  for  acquiaition  and  the  other  dealing  with 
adaptive  technique*  for  eatimating  ayatem  parfcrmance  in  real-time  onbeard  the  vehicle. 


5.  1  Performance  guideline 


A*  pointed  out  above,  tiie  performance  criteria  for  acquiaition  ia  concerned  with  the  avoidance  of  falae 
fUe*  aa  meaaured  by  Ita  probability  of  occurrence,  Pff.  Developing  eom*  general  theoretical  guidalinaa  in 
thia  area  avoid*  the  expenaea  aaaociated  witit  developing  guideline*  completely  from  Monte  Carlo  aimula- 
tiona.  The  gener^  tiieoretical  development  of  determining  Pjf  or  P_  (I  -  Pfj))  begin*  witi>  examining  tite 
correlation  aurface  ahown  in  Figure  6.  The  correlation  value*  can  ^  broken  into  two  graups--thoae  aaeo- 
ciatad  with  match  and  nonmatch  correlation  valuee,  •  (J),  which  are  located  away  from  tiie  central  peak.  A* 
eeen  in  Fig  we  U  theae  correlation  value*  can  be  compactiy  repreaantad  by  two  atatlatical  dlatributiena  - -on* 
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astocUMd  with  th*  nonm&tch  v«lu«i  «nd  on«  asiocUted  «lttt  match  valua(a). *  Tha  correlation  valna(a) 
aaaoeiatad  with  tha  eorract  match  point  may  alao  teka  on  a  dlatribution  of  valuaa  dua  to  noiaa  and  othar 
arrort  (rnch  as  gaomatric)  in  tha  syatam.  Errors  will  hava  a  tandaney  to  spra&d  out  botit  d>a  match  and 
nonmatch  corralation  distribntions.  The  computation  of  tha  parformanca  maasura  involvas  datarmining  tha 
probability  that  a  correlation  valua  drawn  from  tha  distribution  associatad  wittt  tha  match  point  axcaads  all 
corrslation  values  drawn  from  di  .ribution  associated  widi  nonmatch  valuas. 


NON  MATCH  POSITIONS  MATCH  POSITION  NON  MATCH  POSITIONS 


Mathematically  this  can  be  a:>prassad  as  the  general  axvraasion  shown  in  Figure  7,  where  it  is 
necessary  to  compute  for  a  given  match  correlation  value  ttta  probability  tiiat  tha  match  correlation  valua 
axcaede  ^a  nor  r.atch  correlation  value  for  all  nonmatch  pesitione  of  comparison  (Q)  betwaan  the  sanaed 
image  and  raferance  maps  with  this  axprassioa  being  weighted  by  the  distribution  associatad  with  tlie  match 
valuas  (30t  31).  If  tha  match  and  nonmatch  correlation  valuas  are  indeed  indapendent.  dtan,  as  shown  in 
step  2  of  Figure  7,  iha  probability  expression  can  be  computed  using  two  separata  distribution  functions >-000 
for  the  match  value  and  one  for  all  the  nonmatch  val''as. 

In  Ate  real  world  there  are  generally  spatial  patterns  in  the  sensed  imagery  which  are  partially  repeated 
in  some  position  of  the  reference  map.  This  scene  interredundancy  problem  can  be  a  major  source  cd 
system  failures  when  compounded  by  noise  and  otiter  error  sources.  It  also  generally  causes  tits  correla¬ 
tion  value  at  some  nonmatch  points  to  be  highly  dependent  on  the  match  correlation,  dius  preventing  the  two 
distributions  to  bo  separated  and  requiring  a  Joint  distribution  expression  to  be  used  in  computing  the  proba¬ 
bility  that  a  nonmatch  correlation  value  exceeds  tits  match  correlation  value.  If  one  attempts  to  be 
mathematically  correct  in  modeling  this  scene  intorreducdancy  problem,  die  expression  Invol^ng  the  Joint 
distribution  function  (for  match  and  nonmatch  values)  causes  one  into  a  scene  epecific  "modus  opsrandi"  widk 
a  probability  expression  which  iu  too  complicated  to  derive  general  results  f  rom. 


Most  auttkora,  in  attempting  to  develop  a  general  Pc  guideline,  have  ignored  Ate  implications  of  the  scene 
redundancy  and  have  assumed  the  match  and  nonmatch  correlation  values  to  be  independent.  The  implies- 
of  avoiding  to  model  the  Interscene  redundancy  problem  are  twofold.  First,  and  foremost,  the  analysis 


dbna  to  oerreiation  in  the  scene  elements  themselves  several  values  around  Ate  correct  match  peuk  may  be 


which  ;oUowa  to  dotom  '.n*  th«.  Pc  guldoltnoa  ahould  bo  eoaoldorod  o  limitlag  coa*  whora  aolao  aod  odior 
approprUtaiy  modalad  orrora  domlaata  dia  foUura  knodo*  For  aitoatieaa  wnora  tho  aeoao  aaioctlea  pro* 
havo  dooa  o  good  Job  ia  acroaatog  out  tho  aeoao  roduadaaey  foUuro  modo,  dto  oaalyala  could  attU 
provido  uao^  porformonoo  guid  liaoa.  U,  howovor,  ouffteloot  ollort  waa  aot  tnada  la  proporly  aolocHag 
(.•{•ranea  mapa  to  avu$d  aeoao  rodvodaneiaa,  ayatom  porforntaaea  la  llkoly  to  bo  algalfleantly  worao  than 
pradictod  by  thaao  guidollnoa.  Second,  odkor  approxlmatiooa  aod  aaaun^tloaa  beyond  thU  point  take  on  laaa 
algnlflcanto. 


COlWEUTION 
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Roturning  to  Figaro  7,  la  order  to  achiovo  tho  next  level  of  almpllflcation  from  the  general  oxproaalon  it 
la  uoceaaary  to  aaaome  idl  nonmatch  correlation  valuea  are  Independent.  Two  factora  Invalidate  thla 
aaaumption.  First,  the  nature  by  which  aome  computational  algorldima  (1.%.,  Mean  Abaolute  Ditference 
(MAO))  proceaa  dte  map  elementa  leada  to  adjacent  correlatloa  valuea  being  correlated  aad  hence  not 
independent.  Thla  proUem  can  be  overcome  by  die  avoidance  of  thla  type  oi  algorldim.  Another  problem 
ariaea  from  the  fact  that  real  world  acene  elementa  are  almoat  alwaya  correlate  which  leada  to  ^elr 
aaaoclated  correlation  valuea  being  correlated.  Thla  problem  can  be  overcome  by  modeling  die  acene  to  be 
compoaed  of  a  number  of  Independant  elementa  (leaa  dian  tiha  total  numbar  of  acene  element  aatlmating  the 
number  of  equivalent  Independent  elementa  in  die  acene,  aad  uelng  dila  number  In  the  computation  pro- 
ceaa.  Here  not  only  muat  die  acene  be  acaled  by  the  correlation  length  factor,  but  equivalent  acallng  muat 
be  performed  on  the  number  of  nonmatch  comparlaon  polata. 

Further  almpllflcation  of  the  expreaaloa  atep  3  to  atop  4  roqulrea  all  die  noamateh  correlatloa  valuea  to 
be  identically  dlatribotad.  In  general  die  hetorogeneoua  nature  of  acene  atrncture,  Ue. ,  die  acene  being 
compoaed  of  homogeneoua  reglone  with  different  mean  lataaolty  valuea,  can  negate  dila  aaaumption.  The  uae 
of  algorithma  which  tnnd  to  homogeaine  acenea  (auch  aa  the  hybrid  algorldim)  can  overcome  thla  difficulty 
and  make  thla  aaaumption  more  reallatfc. 

Sinco  correlation  valuea  involve  the  aummlng  of  a  large  number  of  random  varlaUea  (aome  combination 
of  dia  acene  elementa)  the  central  limit  dieoretn  can  be  invoked  to  almplify  the  expreaaloa  In  atep  4  farther. 
Thie  aaatnnption  Impliea  dmt  the  dlatrlbutioa  of  die  match  aad  aonmatw  correlatlm  valuea  la  Canaalan. 
rtnaliy,  a  terther  approximatloa  can  be  implied  to  obtain  a  cloeed  form  expreaalon  for  die  performance 


Tc  munmarias,  ona  cannot  davalop  a  uaeful  fanaral  axpraaaioa  toY  computing  Pc  in  tha  proaanca  of 
Bcana  intar radundiancy  prnblamt  and  oihar  axtarnal  arror  aourcaa  auch  at  noiaa  and  diatortiona.  If  one  can 
ignora  Uia  acana  radundancy  pjroblam  by  a  atringant  acana  aalactlon  procaaa,  it  ia  poaaibla,  uaing  a  aariaa 
of  approximationa  and  aaaumptlona,  to  davalop  a  ajcpraaaion  nrhich  yialda  guidalinaa  for  ayatain  par* 
formanca  in  tha  praaonca  of  noxaa  and  othar  arror  arurcaa.  Tha  ramaindar  ^  ttiia  aaction  will  axamina 
davaloping  adaptiva  tachniquaa  for  aatimatlng  ayatam  parformanca  from  tha  corralation  aurfaca  Itaalf. 


In  ordar  to  anaura  miaaion  affactlvanata  and  aafa  warhaad  arming,  a  critaria  to  a^timata  wholhar  tha 
guidanca  updata  ia  aceaptabla  or  not  ia  daalrabla.  Ona  propound  ta^niqua  involvaa  a  voting  logic  widt  thraa 
aucccaaivu  updata  acanaa.  Hara,  &a  datarminad  fix  point  of  two  of  the  thraa  corralatad  acanea  muat  ba 
matched  witiiin  an  aceaptabla  bound;  alae  the  fix  aaquanca  la  rajactad  for  guidanca  updating.  Although 
aimpla  to  implamant  aiM  auitabla  for  uaa  widt  ralativaly  invariant  rafaranca  araaii,  thia  tachniqua  braaba 


down  whan  tha  update  area  ia  mlaaad  altogether;  or  whan  aignificant  varlationB  from  dta  expected  acana 
aignature  axiat  &at  can  not  ba  modeled  b  priorL  Whan  cotmlad  with  the  iaharant  modeling  limitationa  of 
moat  aenaor  operating  banda,  thia  technique  doaa  not  provide  any  indication  of  tha  uncertainty  Of  die 
individual  fixaa. 


'  Two  baaic  tachniquaa  axiat  which  are  capable  of  providing  a  batter  parformac-:  >  \  akimata  of  updata  quality. 
Tha  firat  involvaa  tha  amalyaia  of  tha  diatrlbution  of  dia  raw  aenaor  acana  data.  ’Tndur  conditiona  of  cloud 
cover  or  kurfaca  anow/ica/watar  die  reaulting  Uandard  deviation  of  tha  diatrlbution  v.’lll  often  approach  that 
of  tha  noiaa  equivalent  (apactral)  power  of  tha  aenaor  itaalf,  and  thua  will  ba  aub8tant.aAly  amallar  dian  that 
.1  from  die  unexpected  update  area.  If  a  multimodal  (hiatogram)  diatrlbution  axiata  where  tha  mean  and 
atandard  deviation  of  a  region  are  aubatantially  different  than  expected,  tha  auapact  pointa  can  ba  labeled 
before  further  proceaaing.  When  the  ratio  of  the  number  of  total  imaged  pointa  varaua  thoae  in  tha  auapect 
diatrlbution  region  each  a  predetermined  value,  the  image  can  be  deleted  from  tha  update  or  '/oting  procaaa. 

A  more  reliable  technique  involvaa  computing  tha  correlatioa,  hybrid  or  feature  extraction  aurfaca,  dien 
uaing  propertiea  or  atatiatica  of  that  diatrlbution  for  aatimeting  update  performance.  A  Hat  of  theoe 
techniquea  in  order  of  increaaing  reliability  ia  given  in  Table  11. 

Table  11.  Statlatlcal  Match  Surface  Methoda  for  Eatimatlng  Fix  Quality* 

1.  Main  peak  amplitude 

2.  Main  peak  to  Hr  at  or  higheat  aidelobe 

amplitude  ratio. 

3.  Main  peak  amplitude  va  background  atatiatica  ratio 

4.  Statietica  of  main  peak  va  background 

5.  Above,  with  compeneation  for  intar-point 

acene  correlatioa  and  algorithm  contribution. 

6.  Above,  with  homogeneoua  region  eegmentation 

of  reference  acene. 


eMethoda  in  increaaing  order  of  affectivenaaa 


i 


While  all  alx  approachea  can  be  uaed  with  correlation  or  hybrid  algoridima,  only  tha  firat  three  are 
compatible  widi  feature  matchiew  techniquea.  (Aa  given  in  Table  It,  the  aix^  and  moat  accurate  approach 
for  fix  performance  eatimation  direedy  Inccrporatea  die  hybrid  algorithm. )  An  analytical  relationahip 
between  aurface  atatiatica  and  original  acene  propertiea  may  not  exiat  becauae  of  the  uae  of  edge  or  vertex 
data  for  matching  with  feature  matching  techniquea.  In  either  caae,  a  deciaion  threahold  baaed  on  aurface 
propertiea  or  atatiatica  for  fix  acceptability  miut  be  determined  4  priori  from  accurate  aimulation  of  the 
update  areaa.  The  firat  diree  caaea  are  aimple  to  Implement,  and  utiliae  the  main  peak,  ita  ratio  widi  the 
firat  or  higheat  aidelobe  peak  or  ita  ratio  witt  die  aurface  background  atatiatica. 

The  firat  caae  uaea  the  amplitude  of  the  main  peak  and  ia  generally  unreliable.  It  ia  highly  dependent  on 
die  imaged  "information  content"  (i.  e. ,  Nj)  which  can  vary  aignificantly  with  environment^  diatortiona  and 
SNR.  The  aecoad  utiliaea  die  ratio  of  die  main  to  firat  or  higheat  aidelobe  peaka.  It  ia  generally  unreliable 
unleaa  the  ratio  ia  very  high  or  low.  For  realietic  intermediate  caaea,  die  ratio  will  oacillate  conaiderably 
duo  to  environmental  diatortiona  and  SNR  varlationa.  Although  an  improvement  over  preceding  caaea,  uaing 
die  ratio  of  the  main  peak  veraua  background  atatiatica,  can  ^ten  be  unreliable  becauae  eatimatea  of  the 
original  reference  and  aenaor  acene  atatiatical  propertiea  are  omitted  which  impact  the  matching  aurface, 
hence  diia  ratio. 


The  final  three  approachea  of  varying  dogreee  of  completeneae  determine  a  probability  of  correct  match 

.k a  -  n.^t.-.-ai-ae^a„  naei,  naah  and  hirtiidiii  llumUlil  jMrtgthrtfclM 


Th«  foartti  caie  atillsca  eatlmataB  of  th«  main  paak  and  l>aclcground  atatiatlca  to  datarmina  tha  dagraa  of 
aeparation  bctwaan  in  and  out  of  ragiatar  diatributlona  (banco  P^).  Although  an  improvomant  ovur  tha 
procading  caaaa,  it  doaa  not  utiliaa  aatimataa  of  tha  original  acana  utatiatica  and  ia  aloo  'liaaad  by  ti>a 
matching  algorithm  itaalf..  An  axampla  of  Oiia  caaa  ia  tha  Bhattachryya  diataaca.  In  ttia  naxt  caao,  eom- 
panaation  ia  mada  for  both  ttia  interpoint  acana  correlation  and  tha  impact  of  tiio  matching  algorithma  on  tha 
aurface  ataHatira.  If  level  ahifta  batwean  ragiona  in  ttia  update  area  do  axiat,  thia  approach  can  ha  vary 
accurate. 

Tha  laat  approach  utiliaaa  the  methodology  of  tha  pravioua  caaa,  but  wid»  region  aagmentation  for  hybrid 
proceaaing.  Here,  die  reference  acana  ia  aagmentad  into  homoganaoua  ragiona  and  matched  againat  the 
unaegmanted  aenaor  imago.  Since  diurnal  or  aaaaonal  level  ahifit  variationa  occur  in  moat  apactral  imagery, 
compenaation  to  region  boundariea  ia  generally  roqairad  to  anaura  die  accuracy  of  fix  quality  aatimataa.  The 
hybrid  approach  ia  generally  i-nore  reliable  than  one  which  aagmanta  both  rafaranca  and  aenaor  acanaa  before 
proceaaing;  aince  thia  method  tenda  to  amplify  environmental  and  SNR  induced  region  boundary  diatortiona. 

It  ia  eatimatad  that  thia  hybrid  algorithm  haa  conaidarabla  utility  in  fix  quality  evaluation;  ainca  it  incorpo- 
rataa  tha  atatiatlcal  propertiaa  of  both  tha  original  and  correlated  reference  and  aenaor  acanaa. 

6. 0  Summary 

Thia  paper  advancea  aeven  major  pointa.  Firat,  diara  ara  map  diffaranca  errora  between  tha  rafaranca 
map  and  aenaor  image  which  have  time  and  apatial  varying  componenta.  Tha  magnitude  of  dieaa  errora  ia 
highly  aenaor  wavelength  or  frequency  dependant;  however,  the  atatiatlca  of  die  map  difference  errora  can 
be  quantified  for  each  aeneor  wavelength  aa  a  function  of  the  material  propertiaa  of  the  acana. 

Second,  an  Important  aapect  of  die  problem  it  to  chooaa  and  to  evaluate  reference  mapa  to  avoid  ualag 
areaa  which: 

1)  do  not  contain  aufficient  information, 

2)  have  a  acene  redundant’ ^  problem,  and 

3)  have  materiala  nt  a  wavelength  under  inveatigation  with  large  aignature  oacidatlona. 

Third,  grouping  errora  and  algorithma  into  die  generic  claaaea  indicated  ia  thia  paper  aimplifiaa  tha 
analyaia  and  anablea  the  problem  to  be  atructured. 

Fourth,  certain  algorithma  can  accommodate  certain  claaaea  rf  errora  more  readily  than  odier  typea  of 
algorithma.  Aa  certain  aenaor  wavelengtha  have  a  claaa  of  errora  which  dominate,  it  ia  poaaible  to  pre¬ 
determine  which  algorithm  ia  moat  appropriate  for  dealing  with  acene  data  at  a  given  aenaor  wavelength. 

Thia  paper  definea  thoae  algorithm/aenaor  wavelength  relationahip*  for  aevaral  apecific  operational  condi- 
tiona. 


Fifth,  the  computation  of  the  probability  of  correct  match  ia  acene  dapendenti  and  hence  any  ganaraliaa- 
tion  muat  be  conaidered  an  approximation.  Since  no  abaoluta  meaaurea  can  be  determined,  it  ia  not 
uaeful  to  develop  optimal  algorithma  baaed  on  mathematical  approximationa  to  die  general  P  formulation. 
The  more  appropriate  problem  ie  to  obtain  the  correct  algorithm  for  accommodating  the  map  difference 
errora  which  are  anticipated  to  occur  and  not  to  worry  about  which  aub-claaa  of  algoridim  ia  mathematically 
optimal  for  die  ideal,  nonreal  world  caaa. 

Sixdi,  it  ia  poaaible  to  improve  the  proceaa  of  updating  miaaile  poaition  by  uaing  map-matching  aurface 
data  to  eatimate  in  real-time  the  performance  of  the  eyatam.  Thaae  aatinnataa,  while  approximationa,  have 
proven  tnrough  experimentation  uaeful  in  aeparating  true  matchea  from  falae  matchea  and  can  be  uaed  in 
weighting  the  accuracy  of  the  fix  poaition. 

Seventh,  a  new  claaa  of  map-matching  algorithm,  the  hybrid  algorithm,  waa  preaented  which  incorporatea 
many  of  the  advantagea  aaaociatad  with  the  feature  matching  algorithma  while  avoiding  many  of  die  pitfalla 
aaaociated  widi  extracting  featurea  from  noiay  aenaor  imagea.  It  waa  ahown  to  have  a  aignificant  utility  in 
dealing  with  a  large  number  of  map  difference  errora. 


Concluaiona  and  recommendi 


The  major  atumbling  block  in  analysing  map-matching  ayatema  ia  the  "acene. "  Variationa  in  the 
temporal  and  spatial  characteristics  of  scenes  mitigate  die  need  for  high-o£der  algoridim  refinement  and 
Invalidate  sophisticated  math  modeling  of  the  proceaa.  Such  variationa  ia  scene  imagery  are  the  major 
problems  in  developing  an  automated  aystem.  Two  major  eatitiea  are  required  to  deal  widi  the  problem; 

(1)  the  eatabliahment,  and  (2)  an  analysis  of  a  data  baaa  devoted  exclusively  to  the  image  dynamics  probl  i. 


A  data  base  should  be  created  for  each  aenaor  wavelengdi  enumerated  in  thia  paper.  The  data  base 
should  consist  of  (1)  a  atatiatically  representative  set  of  reference  mapa  covering  Ae  range  of  expected 
materials,  material  interfacea  and  target  types  likely  to  be  encountered,  (2)  aa  accompanying  set  of  sensor 
images  (contained  widiin  the  reference  map  boundary)  which  reflect  the  range  of  expected  temporal  signa¬ 
ture  variations,  and  (3)  a  library’  of  the  physical  and  wavelengdi  dependent  electric^  properties  of  "common" 
acene  materiala. 


Having  davalopvd  auch  •  data  bata  it  la  than  nacaaatry  to  atatlf ti::ally  qiuintlfy  tha  nature  of  the  acene 
errora  preaant.  Baaed  on  ^la  quantification,  the  moat  appropriate  algorithm  claaa  and  preproceaelng 
choice  for  a  particular  wavelength  (and  poaaibly  target  type)  can  be  determined.  Finally,  after  evaluating 
eyctem  performance  over  the  expected  range  of  acene  errora  and  operational  conatralnta  it  would  be  poasi> 
ble  to  determine  which  aeneor  wavelongtha  are  moat  appropriate  for  the  image  matching  taaka. 
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App.'tndix 

Texture 

The  concept  of  texture  from  a  scene  composition /sensor  resolution  description  has  previously  been  intro¬ 
duced.  It  is  of  interest  to  explore  the  source  of  the  textural  spectral  signature  variations  present  on  a  time 
and  space-varying  basis  within  sensor  imagery.  Texture  is  effectively  produced  by  two  different  types  of 
phenomena. 

In  the  first  case,  a  material  of  different  physical  or  electrical  properties  than  its  neighbors  exists  within 
the  scene.  If  it  occupies  one  or  more  scene  elements  it  will  be  resolved  within  the  sensor  image.  It  is  also 
possible  to  resolve  this  material  un  a  subelement  basis  if  its  area  times  radiance  is  greater  than  that  from 
all  other  materials  within  the  element,  and  if  the  resulting  radiance  is  above  the  detectors  SNR. 


In  the  second  case,  a  truly  homogeneous  material  may  exist  within  a  region  virtually  independent  of 
resolution  (i.  e. ,  dry  beach  sand).  Texture  may  still  be  present  due  to  three  factors. 

The  first  is  the  slope  and  slope  azimuth  of  the  material  relative  to  the  sun  (or  iltuminator)-target-sensor 
geometry.  It  is  possible  for  areas  with  moderate  slopes  (20*^— 30*^)  to  produce  substantially  more  or  less 
radiance  depending  on  their  orientation  between  illumination  source  and  detector  (this  is  especially  impor¬ 
tant  when  direct  (i.  e. ,  laser  or  sun)  versus  diffuse  (i.  e. ,  skylight)  irradiance  exists). 

For  imaging  lasers  and  radars,  slope  and  slope  azimuth  between  the  illuminator-target-detector 
geometry  significantly  impacts  the  magnitude  of  the  returned  target  radiance.  At  least  here,  versus  a 
passive  system,  the  illuminator  is  usually  co-located  near  the  detector.  The  net  effect  of  this  is  to  simplify 
the  governing  geometry  for  determining  the  return  vector  of  tlie  propogated  wave.  Weak  to  moderate 
reflectors  oriented  at  steep  incident  angles  to  a  co-located  transmitter /detector  can  often  produce  a  signifi¬ 
cantly  greater  return  than  strong  reflectors  oriented  less  favorably. 

For  passive  systems  where  the  illuminator  (usually  the  sun)-vehicle  geometry  is  generally  not  co-located, 
shadowing  is  more  difficult  to  evaluate.  Here,  shadowing  is  often  a  problem  due  to  the  time-varying  sun- 
target  geometry  coppled  with  slope  and  slope  azimuth,  aud  surface  roughness  of  the  reference  area.  As  a 
consequence,  shadowing  can  significantly  impact  daytime  optical/near  IR  and  middle  IR  imagery  where  a 
strong  solar  component  exists.  Its  effect  in  the  thermal  IR  region  is  to  prevent  direct  incident  short  wave¬ 
length  radiation  from  being  absorbed  by  the  target,  thus  reducing  the  diurnal  temperature  oscillation  by 
weakening  the  thermal  inertia  driving  function.  Because  of  the  diffuse  nature  of  passive  microwave  radi¬ 
ation,  the  effect  of  shadowing  on  apparent  brightness  temperature  is  generally  not  a  problem  within  a 
reasonable  range  of  antenna  depression  angles  with  this  form  of  imagery  for  water  and  metal  because  of 
their  moderate  to  high  microwave  reflectances  respectively.  For  materials  with  high  microwave  emittances, 
shadowing  weakens  the  thermal  inertia  driving  funct'on  for  diurnal  temperature  oscillation,  and  can  reduce 
the  emissive  power  and  the  resulting  observed  apparent  brightness  temperature. 

As  with  active  illuminator  systems,  slope  and  slope  azimuth  play  an  important  role  in  many  passive 
imaging  systems.  In  the  optical/near  IR  and  middle  IR  regions  it  impacts  the  returned  target  radiance 
similarly  to  active  systems, although  to  a  greater  extent  because  of  the  varying  solar-target  geometry.  Slope 
and  slope  azimuth  also  produce  different!^  heating  from  absorbed  short  wavelength  solar  radiation,  which 
can  have  a  moderate  to  strong  impact  on  night-time  middle  IR  and  diurnal  thermal  IR  imagery  and  a  small  to 
moderate  effect  on  diurnal  passive  microwave  imagery  when  a  high  microwave  material  emittance  exists. 

The  necond  parameter  that  can  produce  image  texture  is  the  material  reflection  coefficient  or  reflectanca 
A  variati'^n  in  smooth  surface  reflected  energy  versus  incid:.nce  angle  exists  due  to  (real  and  imaginary) 
material  elec‘'rical  components.  Because  material  reflectance  is  the  dominant  energy  balance  parameter  in 
several  spectrv'  regions,  its  directional  characteristics  can  have  a  significant  impact  on  the  amount  of 
energy  returned  irom  a  target.  The  real  component  is  the  material  dielectric  component,  while  the  imagi¬ 
nary  one  equals  th*>  electrical  conductivity  divided  by  the  angular  frequency  times  the  free-space 
permittivity.  For  conductors  (i.  e. ,  bare  metals),  tire  second  component  generally  predominates.  For 
dielectrics  (most  nature!  materials)  the  first  term  is  usually  dominant  since  negligible  electrical  conduc¬ 
tivity  exists. 

Given  the  electrical  component  values,  the  vertically  or  horizontally  polarized  directional  reflectance  for 
a  smooth  material  can  be  determined  from  Fresnel's  equations  at  a  particular  wavelength  (34).  The  values 
computed  by  Fresnel's  equations  provide  a  measure  of  theoretical  material  reflectance  versus  incident  and 
reflected  ‘  'gles.  Surface  roughness  height  and  orientation  can,  however,  significantly  impact  the  amount  of 
radianr  tually  reflected  (or  absorbed)  by  the  material. 

Consequently,  the  third  paramwter  of  interest  is  the  roughness  of  the  surface  itself.  For  the  same 
illumination  source-detector  geometry  multiple  reflections  will  occur  within  the  material  when  the  roughness 
height  to  wavelength  ratio  is  large.  This  results  in  an  increase  in  absorptance  and  a  consequent  increase  in 
emittance  in  wavelength  regions  where  this  parameter  is  relevant.  When  the  ratio  of  roughness  height  to 
wavelength  is  small,  multiple  reflection  effects  diminish,  and  the  resulting  absorptance  decreases  to  a 
theoretical  minimum  for  the  material. 
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The  directional  and  bidirectional  reflectancea  of  a  material  in  a  given  wavelength  region  are  dependent  on 
the  surface  roughneaa  a*  well  as  the  governing  electrical  relationship.  Unfortunately,  it  is  difficiQt  to 
characterise  surface  roughness.  The  root  mean  square  (rms)  roughness  sometimes  used  provides  no  infor- 
mation  pertaining  to  the  statistical  distribution  of  roughness  around  the  rms  value  and  the  average  slope  of 
the  sides  of  the  roughness  peaks,  although  tiiey  can  significantly  influence  directional  and  bidirectional 
reflectance  (34).  As  a  result  the  directional  and  bidirectional  reflectance  characteristics  dtat  some 
materials  ej^iblt  are  due  to  a  combination  of  complex  surface  structure  and  surface  impurities  (i.  e. ,  iron 
oxide  or  moss)  present  together  with  the  materials  inherent  electrical  properties. 

Surface  roughness  significantly  impacts  the  returns  from  Imaging  lasers  and  radars  depending  on  the 
relevant  geometry.  When  a  shallow  incidence  angle  exists  between  the  illuminator  and  target,  the  net 
effect  for  rough  surfaces  is  often  to  return  more  radiance  than  from  a  smooth  one  because  of  the  effective 
presence  of  material  corner  reflectors.  This  is  evident  in  radar  data  when  examining  imagery  from  smooth 
versus  rough  fields  or  water.  At  steep  incidence  angles,  flie  reverse  is  true.  Here,  a  smooth,  surface  will 
generally  return  more  radiance  than  from  a  rough  one. 

Surface  roughness  can  also  impact  the  returned  target  radiance  in  the  optical/near  IR,  and  to  a  lesser 
extent  in  the  middle  IR  region,  because  of  the  predominance  of  the  direct  solar  illumination  component.  The 
effects  are  similar  to  those  discussed  under  active  illuminator  systems.  In  the  thermal  IR  region,  an  insig¬ 
nificant  amount  of  direct  energy  radiated  by  the  sun  reaches  the  surface  and  generally  high  emittance  exist 
for  most  natural  materials  (typically  .  85  to  .99).  As  a  consequence,  the  net  effect  here  is  to  impact  the 
short  wavelength  absorptance  and  possibly  the  material  thermal  Uiertia. 

In  the  passive  microwave  imagery  where  high  emittance  materials  are  present  which  have  a  surface 
roughness  substantially  greater  than  the  Imaging  wavelength,  effects  similar  to  those  in  the  thermal  IR 
can  exist.  Many  materials  such  as  metal,  concrete,  asphalt  and  smooflt  water  behave  specularly  in  the 
passive  microwave  region  (particularly  at  frequencies  below\40GHz)  because  deir  surface  roughness  is 
small  in  comparison  'jo  the  imaging  wavelength. '  An  interesting  case  of  the  effect  of  surface  roughness  on 
material  reflectance  in  this  region  occurs  with  water.  Calm  water  behaves  as  a  good  specular  reflector  of 
passive  microwave  radiation  (second  only  to  metal  in  this  wavelength  region).  As  surface  roughness 
increases,  the  magnitude  of  tlie  sky  radiation  times  microwave  material  reflectance  term  decreases  due  to 
multiple  reflections  present.  As  a  consequence,  die  emittance  times  the  ground  (water)  temperature  term 
predominates  in  rough  water  where  capillary  waves  exist,  and  emissive  power  variations  under  clear  skies 
can  be  on  the  order  of  20ft  to  SOfc  between  this  and  the  smooth  water  case. 


